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DYNAMICAL EFFECTS OF X-RAY INTERFERENCE IN NATURAL 
NaCl CRYSTALS FROM WIELICZKA 


By J, AULEYTNER AND M. PAWŁOWSKA 
Institute of Experimental Physics, Warsaw University 


(Received January 2, 1961) 


The dynamical theory of X-ray diffraction is based on a consideration of the interaction 
between the primary and diffracted waves. This theory was developed by Ewald (1933), 
Laue (1949, 1952), and Kato (1952). 

Borrmann (1955) showed experimentally that interference affects the intensity of the 
primary beam propagating through the erystal. Further investigation of these effects brought 
to light the effect of anomalous absorption of X-rays which, on being incident on the crystal 
under the angle fulfilling Bragg's condition, are diffracted and traverse the entire bulk of 
the crystal. This discovery opened up a series of experimental and theoretical investigations 
dealing primarily with quasi-ideal crystals, (Brogren, Adell 1954), (Guinier, Guyon 1959). 

The interaction between the primary and diffracted waves within the framework of the 
dynamical theory can be represented by internal reflections within the crystal, as shown 
schematically in Fig. 1. 

In this respect, it is necessary to distinguish between a) weakly absorbing crystals and 
b) strongly absorbing crystals. 

In crystals that are weak absorbers of X-rays, the so-called internal reflections occur 
within the entire region admissible for a given interference. As the beams subject to reflec- 
tion are weakly absorbed, they can give rise to interference throughout the entire region. 
For strongly absorbing crystals, these internal interactions are shown schematically in 
Fig. lb. The reflections can be said to take place within a limited region, as the beam is 
strongly absorbed. Only pencils for which the distances between the consecutive reflections 
are small reach the surface and can emerge from the crystal. 

A more exact analysis of the manner of propagation of the X-ray beam in the crystal 
involves the problem of the energy flow therein. The theory based on the laws of propaga- 
tion of an electromagnetic wave in a medium enables us to compute the mean value of the 
Poynting vector (the vector of energy flow in the wave field). 

According to the dynamical theory, if Bragg's condition is fulfilled, absorption within 
the crystal depends on the direction; it is the smallest in the direction of the difiracting 
planes. Fig. 2 shows the transmission curve for beams subject to interference within the crys- 
tal, being the locus of the points reached by the energy of the incident beam on attenuation 
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Fig. 1. Schematical representation of interactions between primary and diffracted waves: 


a) for weakly absorbing crystals, b) for strongly absorbing crystals. 


in the proportion of e 1. It is seen that the least degree of weakening occurs in the case of 
pencils propagating in the direction of the reflecting lattice planes. 


The foregoing cases of strongly absorbing crystals were experimentally investigated 
by Borrmann and Guinier. 


Borrmann experimented on calcite crystals strongly absorbing the characteristic molyb- 
denum radiation applied. The region of interference within the crystal was very narrow, so 
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that the following four reflections could be distinguished: the first, due to the transmitted 
beam, the second having the direction of normal Bragg diffraction, the third arising from 
internal interference and emerging in the direction of the primary beam, and the fourth, 
also due to internal interference, in the direction of the diffracted pencils. 

Guinier investigated mosaic and weakly absorbing crystals. On the photographic plate, 
in the immediate vicinity of the transmitted beam, he observed a blackening due to internal 
reflections of the X-ray beam within the crystal. It had been Guinier's assumption that the 


Fig. 2. 


crystal is not required to be an ideal one in order that the effect of dynamical scattering 
should appear in weakly absorbing crystals. For this, it is sufficient that small regions distrib- 
uted throughout the entire bulk of the crystal, having identical lattice planes, mutually 
parallel, and fulfilling Bragg's equation in given conditions, shall be present. The reflections 
obtained by Guinier exhibited a discontinuous structure. This was due to defects in the 
crystal. 

It was the aim of the present investigation to verify whether dynamical X-ray effects 
can occur in the case of NaCl. The hypothesis arose that certain specimens of NaCl crystals 
from Wieliczka may present regions of quasi-ideal structure, thus providing the case of a crys- 
tal moderately absorbing the characteristic Mo radiation. 

The investigation was carried out with the system shown in Fig. 3. The source con- 
sisted of an X-ray bulb with molybdenum plate. The characteristic radiation as obtained by 
interposing zirconium filters in the path of the incident beam was utilized. On passing through 
the collimator, the beam reached the crystal, which was set up on the spectrometer support. 
The crystal was placed in such a position that the X-rays were reflected from planes perpen- 
dicular to the crystal surface (in the present case, from the planes (100)). The photographic 
method was used for detecting the effect. The photographic plate was placed behind the 
crystal, parallel to the latter. The pirmary beam was screened off with a Pb plate in order 
to make possible the detection of any X-ray scattering in the direction of the primary beam. 
The Pb plate was interposed between the crystal and the photographic plate, parallel to the 
crystal. The Geiger-Müller counter L was used for verifying whether the the position of the 
crystal was such as to yield the maximum of the interference fringe (the G-M counter II served 
for testing stability and measuring the intensity of the beam emerging from the bulb during 
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measurements of the coefficient of absorption). The majority of the crystals investigated 
possessed a structure sufficiently ideal for the observation of dynamical effects. Fig. 4 
the picture of the effects as obtained. with reflection of the X-ray beam from the family 
of planes (100). To the left, we see a weak line, with the blackening inhomogeneously distribą 
uted, due to dynamical scattering of the X-ray beam in the direction of the transmitted 
pencil. The transmitted pencil itself had been screened off (the screen being to the left of 
the line under consideration). The more intense, bright light band to the right arose through 
reflection of the incident beam from the family of planes (100). The left hand portion of 
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Fig. 3. 


this band is somewhat lighter. This is the portion due to dynamical scattering in the direction 
of the reflected beam. The central, oblique fringe is due to reflection from the planes (111). 
The reason for its presence resides in a degree of divergence in the incident beam. Fig. 5 
shows the trace of the X-ray beam transmitted by a crystal placed in such a manner that 
the only reflection should be from the planes (111) (the more intense band on the right 
hand side of the picture). In this position of the erystal, we have also the effect of dynamical 
scattering in the direction of the transmitted beam (the weak, bright line to the left) and 
in that of the reflected one. This picture shows even more clearly the difference in intensity 
between the portion of the reflected beam related to the dynamical scattering and the re- 
flected beam proper (the right hand side of the Bragg fringe is weaker). However, the band 
is not divided, as was the case in Borrmann’s experiments on calcite crystals. The intense 
bright fringe on the left hand side of the weaker band is the trace of that portion of the in- 
cident beam that was not screened off. The fact that dynamical scattering occurs also on the 
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Fig. 6. Photograph of NaCl crystal surface taken on etching in aqueous alcohol solution 
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planes (111) when these are the only ones to fulfill Bragg’s condition suggested that scattering 
in the direction of the incident beam may represent a superposition of scattering on the 
_ planes (100) and (111), if reflection occurs simultaneously from either family of planes. 
. Guinier computed the intensity to be expected in a beam resulting from dynamical scattering 
in the direction of the transmitted beam, for the case of interaction with a single family of 
planes. The theoretical result is in good agreement with the intensity measured. Scattering 
in the direction of the transmitted beam would seem to be a superposition of scattering by 
the families of planes (100) and (111), if there is simultaneous reflection from either of these. 

Etching tests proved the natural crystals already mentioned to possess a quasi-ideal 
structure. Fig. 6 is a microscopical picture of the surface of an NaCl crystal subsequent to 
etching in an aqueous alcohol solution. The dislocation density amounts here to approxi- 
mately 300 cm-?. 

Thermal etching also yielded a dislocation density of the same order. Heating the spe- 
cimen at about 750? C for an hour in vacuum of the order of 1074 mm Hg failed to produce 
any changes in the X-ray diffraction pictures. 

The X-ray interference photographs obtained reveal, to a certain degree, the inner 
structure of the crystals. Irregularities in the lattice are revealed by the sharp rise produced 
in the absorption of the interfering beams (the reflections shown in Figs. 4 and 5 are 
not of equal blackening). 

From a comparison of the present authors' results with those of Guinier, the natural 
NaCl crystals are seen to yield a pattern similar to those obtained in the case of weakly 
absorbing crystals. 

Investigation of the intensity distribution in the reflections will continue with the aim 
of checking with the theoretical results of Laue for NaCl. Investigation of the intensity re- 
partition requires a strictly parallel beam, in order to avoid superposition of the effects 
of dynamical scattering on reflection from two lattice planes. 


Conclusions 

1. Dynamical effects of X-ray interference were found to occur in natural NaCl single 
crystals from Wieliczka in the case of reflections (200) and (111). 

2. The crystals exhibiting these effects reveal a relatively low dislocation density, which 
explains the fact that the effects were accessible to observation. 

3. Heating the specimes failed to produce changes in the X-ray diffraction patterns, 
thus pointing to high stability of the structure. 

The authors are greatly indebted to Professor L. Sosnowski for directing the present 
investigation and for his discussion of the results. 
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YANG Kuo-SHEN 
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Using the self-consistency condition and the first order perturbation method and assuming that the nuclear 
self-consistent field is the Nilsson potential, the deformability of nuclear cores and the electric quadrupole moment . 
of nuclei with a doubly closed shell core + one nucleon are calculated. The former case shows exact agreement 
with observation and the latter fits the experimental data as well as other best methods. 


1. Introduction 


The large electric quadrupole moment and the relatively low half-life of nuclei with 
a doubly closed shell core + one nucleon, e. g. 01", F", Pb?®’, cannot be explained by the 
single particle model with spherically symmetric potential. The electric quadrupole moment 
caleulated by the nucleon-surface coupling model fits the experimental data rather well 
for some nuclei (Barker 1956, True 1956, Raz 1957, Mukherjee and Dutt 1958). However, 
in this calculation the observed half-life and energy difference of the E2 transition have 
been utilized. Theoretically this model does not take into account the antisymmetry require- 
ment between the external nucleon, and the nucleons of the core, overestimates the con- 
tribution from the surface nucleons and underestimates that from the interior nucleons. 
Some authors attacked the problem by means of the method of interconfigurational mixing 
(Perks 1955, Horie and Arima 1955, Amado and Blin-Stoyle 1957, Barton 1959, Barton, 
Brick and Delves 1959/60); some others, used both the nucleon-surface coupling model 
and the interconfigurational mixing method (True and Ford 1958). They got considerable 
success for a very small number of nuclei, although some difficulties exist in their theoretical 


bases (Barton 1959, Barton, Brink and Delves 1959/60). Besides some authors have estimated 
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with correct order of magnitude! the electric quadrupole moment of 077 on the basis of 


Brueckner's theory (e. g. Choong 1959). E 
Another method (Mottelson 1958) is based on the assumption of deformability of the 


nuclear core by the external nucleon or hole. The total quadrupole moment Qy, of the nucleus 
is the sum of two parts: 

On = o = Q5 
with +(—) sign for an external nucleon (hole). The deformability of the nuclear core is 


defined by the relation: 


Geom re QY (1) 
The intrinsic electric quadrupole moment may be written as 
Kaa. eZ. £ 
Qi ey Q, FE Qoy 2) 
Pong 


where A, and Z, denote respectively the mass number and the number of protons of the 
core: ep=0 for an external neutron or neutron hole. and e,=e for an external proton or 
proton hole. 

Experimentally, y= 1.0 for OY and F”, and y=2.9 for heavy nuclei in the neighbour- 
hood of Pb? (True 1956, Szymansky 1959). Theoretically, y may be calculated by means 
of the first order perturbation method on assuming that the nucleus deforms self-consistently, 
i. e. the nuclear shape determined by the wave function coincides with the shape of the 
potential well. An appropriate choice of the potential well is very important. For harmonic 
oscillator potentials the calculated y= 0.82) for OY and F”, and y =1 for heavy nuclei (Mottel- 
son 1958). For a square well of infinite depth, the calculated y>3.5 for all nuclei (Szymansky 
1959). In this paper a much more realistic potential, the Nilsson potential is chosen. The 
assumption that the nucleus deforms self-consistently will be used to determine not only 
the deformability but also the equilibrium deformation and the electric quadrupole moment 
of nuclei with a core + one nucleon, 


2. Methods and Formulae 
Suppose the nuclear shape is aa ellipsoid of revolution 
r = Ro(1+ dP,(cos 0"), 


where r”, ©”, ®’ and x,’ y', z' are the spherical and Cartesian coordinates respectively in the 
body fixed frame, and d is the deformation parameter related to 6 (Nilsson 1955) and p (Bohr 
and Mottelson 1953) by the relation 


a= = — p 


! The estimated order of magnitude is —]1 x 10% em? for 017, 


f * For light nuclei the value y=1 reported by other authors, is obtained by neglecting, in the calculation 
of b (see Eq. (10) below), the contribution from the external nucleon or hole. The authors wish to thank Dr Z. 
Szymansky for pointing out this neglection in a kind letter to us. í 
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The condition that the nucleus deforms self-consistently may be put into the form (to first 
order in d): 


SOIN eC: =D) + < pplz” — xz | p) 
"| QNM) s 


where © and g, are the intrinsic wave functions of the nuclear core and the external nucleon 
or hole respectively, >; denotes summation over all nucleons of the core, and --(—) sign 


1 
in the symbol “ +” is taken for an external nucleon (hole). It is noted that Eq. (3) has been 
obtained by means of the principle of minimum energy (Moszkowsky 1956). 
Suppose the self-consistent field is the Nilsson potential (Nilsson 1955). Using Nilsson’s 
notation the single particle Hamiltonian may be written as 


H =H,+ H, 
Hy = &hoq(— 4 + y?) + CL: S+ DB, 
H, = — dhogr?Pą(cos 0). (4) 
The eigenfunctions of H, are denoted by |NZQ> with the eigenvalues, 
Eny = (N+ $)hox — [+ 1) — +1) — š+ C+ Dulha. (5) 


Applying the many-body perturbation method on the nuclear core, © may be expanded 
as a power series in d. All wave functions of the unperturbed system are Slater determinants 
consisting of |NIjQ). A simplified form of the coefficients of the terms belonging to the 
first power of d is: 
hoy CNPJ O |r2P, (cos 0) ENGR 

Eny — Enry 


(6) , 


Cn17’2 


where |N/jQ) is the neutron (proton) state belonging to the core, and |N’/’7’Q’) is the neutron 
(proton) state not belonging to the core. In general d is small for nuclei with a doubly closed 
shell core tone nucleon. If all coefficients Cyr; o are small, the first order perturbation meth- 
od is sufficient. From the numerator of Eq. (6), we see that Cyro #0 only if N' = N or N 
=N +2. Evidently Cyg is small for N' — N +2. If both neutrons of the core occupy all states 
with N< some integer N, and the protons of the core occupy those with N < some integer Na, 
there exists no coefficient with N'=N and the first order perturbation method is sufficient. 
Nuclear cores satisfying this condition are „Het, 301%, „Ca? and ,,Pb?98, Whenever coupling 
coefficients exist with N’=N, the first order perturbation method may still be approximately 
applied for our purposes, if only a small number of coefficients Cyg are large and their 
total contribution to the first power of d is not large. All nuclei of this type are collected in 
Table I. For other nuclei, where this first order perturbation method breaks down, another 
method will be introduced at the end of this section. 
Substituting the above first-order expansion of ®, we get after some simplifications 

= 25) (5° — a.) =2ad (7) 

core 

Heko (8) 

b—a 


MUR 
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where 


a= TY = ADA WICZCOCOO 0 


Mol Eny = Ewry)] h Mo 
Nij N'lj 


h 3 
b= ) «NIjO|r?| NO» = —— ) Nha- (10) 

Mog 2 

core--1 core+1 
X, denotes summation over all, possible values of N, l, j of neutrons and protons belonging 
Nij : ° 7h yen et 
to the core, and >)’ denotes summation over all possible values of N’, l’, j’ of neutrons (protons) 
N 

not belonging to the core whenever N, /, j denote quantum numbers of neutrons (protons). 
The value of h| Me is determined by the observed equivalent radius of the proton distri- 


% 
bution, i. e. «7| = 1.20 x 10-13 45 cm (Hill 1957). Thus we get, 


zZ 
E 
u z 5 20x10 az] M (w+ 3), (11) 
0 "RY 


i=1 


Substituting Eqs. (7), (8) in (1), we get 


a 


ant es a 


(12) 


Eq. (12) tells us that y may be calculated even though we do not know both the equilibrium 
deformation d and the value of h/May. 

The equilibrium deformation d of the self-consistent field can be immediately solved 
from eq. (8), if the explicit dependence of Q, on d is obtained by applying the first order 
perturbation method on the single-particle wave function of the external nucleon or hole. 
For a great majority of nuclei witha core + one nucleon, the first order perturbation method 
with |V/jQ) as unperturbed wave functions can be applied sufficiently well, since in the 
coefficients of d, no coupling terms between states with the same N can appear due to the 
fact that for these nuclei the quantum number Q of the external nucleon or hole is the largest 
among all states with the same N. For other nuclei, where this condition is not fulfilled, the 
modified first order perturbation method introduced at the end of this section will be used. 

To calculate the equilibrium deformation of the ground state of nuclei with a core 4- 
one nucleon, it is necessary to know the magnitude of Q of the external nucleon or hole. 
Theoretically it may be obtained by the principle of least energy. In this paper it will be 
determined by the empirical rule: Img. 

The calculated intrinsic electric quadrupole moment Qa i$ obtained by substituting 
the relation Q, = 2da/y in Eq. (2) and using Eqs. (9)—(12). For nuclei of weak deformation, 
on observing some formula in § V. b) of the paper of Bohr and Mottelson (1953), one can see 
that the intrinsic electric quadrupole moment is approximately equal to the electric quadru- 
pole moment. Thus the calculated intrinsic electric quadrupole moment will be used to 
compare with the observed electric quadrupole moment Oi: | 
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For a harmonic oscillator potential, according to the Thomas-Fermi model, the observed 
deformation parameter dps or 85,, of the proton distribution may by derived from the ob- 
served electric quadrupole moment Q, by the relation (Tseng et al. 1959): 


_100Q 320 
out obs obs 
obs 2Z'ls | Ë = Z 


) (weak deformation), (13) 


The calculated deformation parameter f? ie may be obtained from its relation with d of 
the self-consistent field. Thus we get 


(14) 


As mentioned above, the first order pertubation method with |N 1j@> as unperturbed 
single particle wave functions breaks down, whenever the contribution from the coefficients 
that couple the neutron (or proton) states with a common M, is large. In this case the first 
order perturbation method will be modified by an appropriate choice of the unperturbed 
single particle wave functions, so that all coefficients coupling the neutron (or proton) states 
with a common M, vanish. It is well-known that the Nilsson wave functions corresponding 
to the deformation d possess this property. The difference between the Nilsson Hamiltonian 
and the.unperturbed Hamiltonian with Nilsson wave functions as its exact eigenfunctions 
will be used as the new perturbation. Expanding © by the correspondingly modified many- 
-body perturbation method, we get coefficients similar to eq. (6) in the first order terms. 
For N’ = N, these coefficients vanish. For N’# N, these coefficients are equal to eq. (6) to ' 
first order of d, since as d — 0 the Nilsson wave functions, eigenvalues, and the matrix 
elements of the new perturbation reduce to |NIjQ), Eny and the matrix elements of 
H; respectively. 

Formulae (7), (8) and (11) with |N} Q> simply replaced by the Nilsson wave functions, 
do not hold true in general?. Nevertheless, d may be found directly from Eq. (3) and by 
substituting it in Eqs. (1) (2), y and Qar may be obtained. As an example, we shall apply 
this method to the nucleus ;B. 


3. Results 


The factor >; IL 7 2| Ps(cos 0)|j2>|? in Eq. (9) are most generally evaluated by utilizing 
8 


the Wigner 3j and 6j symbols (see Appendix I). Another factor <N'lr?|NI> of Eq. (9) was 
evaluated by Nilsson (1955). Eq. (5) is used to obtain the factor (Eny — Enq,)[hog, where 
u, X are selected according to the opinion of some authors (Nilsson 1955, Newton 1957, 


Bishop 1960, eżc.): 


3 For the special cores, ,O19, ‚Het, ‚Ca? and ,,Pb?98, this formula still hold true. 


(neutron) (proton) 
» 0, I 2 3 4 | 1 5,6 7 | 8 
X 0.125 | 0.125 0.05 0.05 0.05 0.05 | 0.05 0.05 
n 0 | 0.20 0.35 0.45 0.55 0.45 1 0.40 0.35 


For N < 2, relatively large errors of y and u may arise. However, their effect on the calculated 
Yan: Beate and Qu, is very small. For example, when the core just occupies all states of 
N=0 or those of N=0 and N=1, the calculated values of y, f, and Q, remain unchanged 
by varying x from 0.05 to 0.125 and keeping u=0. For OY, where the core occupies all states 
with N=0, 1 and 2, the calculated values of y are 1.00, 1.02 and 1.04 for 5 —0.08. u=0.17, 
470.125, u=0.20, and 4 —0.125, u=0.33, respectively. For N—8, the selected value of 
u may possess some error. However, the effect of this error to reasonable extent on the calcu- 
lated results is not large, since the calculated y for Pb29 are 3.07 and 2.89 corresponding to, 
for N=8, x=0.05, u=0.40 and y=0.05, u=0.35 respectively. 


TABLE 1 
The deformability, the electric quadrupole moment and the equilibrium deformation 


Nucleus | Yealc. Yobs. Qcale. Qobs. | Beate. | D 

Had 201 | 

„Hej 2.81 

„He3 0.64. 

s Lis 0.64 

„Ni? 1.47 

Ko 1.47 

Oo" 1.02 1.0 —0.027 -0.0265-1-0.003 — 0.046 —0.041 
A 1.02 1.0 

19918 5 97 0.058 0.054 +0.01 0.030 0.026 
1891 2.27 

17C156 2.34 -0.095 | -0.063 0.002 —0.052 —0.028 
19 K30 1.47 0.080 0.09 +0.02 0.034 0.033 
mi 1.47 | 

ao Cal 1.26 

119640 1.26 

"ECC 3.64 

102119 3.69 

soZ0gi 3.22 —0.22 —0.46 —0.023 —0.048 
41Nb50 3.22 

ern 2.50 

505071 2.31 

AART 2.32 —0.37 —0.53 +0.10 —0.025 —0.037 
TIME 2.89 

aPbiez 2.89 2.9 

mura 2.15 

Bilde | 2.76 —0.61 —04 —0.019 —0.013 
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All nuclei of Table 1 are calculated by means of the methods and formule of $ 2. The 
electric quadrupole moment and the equilibrium deformation of all nuclei in Table 1 can 
be easily calculated, although only those of a few nuclei for which we possess experimental 

data, are tabulated. This makes the present method more convenient in application than 
other methods. Comparing Year with the observed value, we conclude that the agreement 
is exact. Since the observed y are derived from the observed half-lives, it appears natural 
that the relatively large half-lives of nuclei such as 07, F", Pb207 etc. might be well explained 
along the theoretical lines of the present method. In general the calculated Qa, and ff. 
fit the experimental data as well as those obtained by other best methods, such as the nucleon- 
-surface coupling model, the interconfigurational mixing method, etc. The larger discrepancy 
for the nuclei 44Z:2; and z,Sbzg', is due to the worse applicability of the first order 
perturbation method. 

Lastly, consider the nucleus ;B!!, Applying the same first order perturbation method as that 
for the nuclei in Table 1, we get b <a, resulting in a large negative value of y. As mentioned in 
82 this is due to the breaking down of this method. Atelrnatively, applying the modified 
first order perturbation method as used at the end of § 2, we get Quae = 0.052, BP. = 0.18 
and Ya. 95-1, to be compared with the observed values: O ps = 0.036, 6%,,= 0.12. The 

discrepancy could be reduced by taking into account the residual force. Nevertheless, this 
discrepancy is of the same order of magnitude as obtained by the nucleon-surface coupling 
model (Mukherjee and Dutt 1958), or by the intermediate coupling model. 


APPENDIX I 
The evaluation of the expression 
Xj |<ety'olp, cos 9 lO: 
Introducing the tensor operator of 2nd rank T(2), the q = 0 component of which is 
T(2, 0) — P, (cos 0), we get by means of some well-known formulae (see Edmonds 1957): 


urere oia = cw (Fo À) WIEN. 


ame = c- ys eren GEY Cro), 


emen - c (o) Teom. 


With the of aid these relations and of the orthogonality properties of the 3j symbol, we get. 


X [X17 Q|P,(cos 0) |j |? = $ G^. j, 2) TOP 


Ki /1 
WU 


? | ("1T (2,0)]10) |? 


KALBA 
—110 


AB 


5 Oj’, J, 2) (A+) (2) + 1) 
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The evaluation of the last expression gives Table 2: 


TABLE 2 
ów e DKR P, (cos 0)|UQ | 
Q 
I, I+ 1/2 1-12 
m Pu 
j = í 
"e 2 Ul+1) (+2) 6 i+ 1) 
ig 5 (21+ 1) (21+ 3) 5 (21—1) (21 + 1) (20 4- 3) 
Be 6 KI--1) 2 (-01(4- 1) 
d BE 5 (21—1) (2141) (2143) 5 (21—1) (213-1) 
For P—14-2 
: RA 


er 1+ 1/2 Rt. 
n i 
1 


3.01) (0452) (13503) 
5 (204-3) (21-45) 


E 6 (I+ 1) (1-2) 3 14-1) (04-2) 
T 5 (21+ 1) (21+ 3) (21 4- 5) 5 (2014-1) QI 4- 3) 
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ON A NEW METHOD OF DETERMINING THE FOCAL PLANE OF 
A SPECTROGRAPH 


By B. Nowak 
Physics Department, Nicholas Copernicus University, Toruń 
(Received April 21, 1961) 


A method for the exact determination of the focal plane of a spectrograph is proposed, utilizing 
a special diaphragm with six slits situated behind the collimator lens. This diaphragm splits the 
image of each defocused spectral line into two groups of narrow interference fringes presenting 
a characteristic distribution of the central intensities enabling to distinguish, within each 
group, the central "marked" fringes. The focal plane is determined and the plane of the plate 
placed thereto by measuring the distance between the "marked" fringes obtained at different values 
of the variable parameters of focusing (the shift of the slit, that of the camera and the angle of 
rotation of the plate holder). Cases of various numbers of focusing parameters (l, 2, 3) are 
discussed and an example is given of the focusing of a /SP-67 spectrograph by means of the dia- 
phragm mentioned. The method of choosing the optimal parameters of the diaphragm for a given 
spectrograph is described in detail. In the Appendix, formulae for the intensity distribution in 
different planes of observation, for the case of Fraunhofer diffraction on diaphragms of the type 
described, are derived and discussed. These are then applied to computing the intensity distribution 
obtained with this diaphragm in several planes. The theoretical distribution is compared with 
the experimental results. 


Introduction 


The utility of a spectrograph with several exchangeable cameras for obtaining spectra 
of different degrees of dispersion is strongly dependent on the possibility of re-focusing the 
spectrograph rapidly after a new camera has been introduced. In general, focusing of a spec- 
trograph consists in finding the values of three parameters, namely, 1 — the shift of the 
slit, L — that of the camera and p — the rotation angle of the plate holder’. The methods 
discussed in the textbooks of spectroscopy (see, e.g. Prokof'ev, 1951, p. 183) consist essen- 
tially in making a number of photographs of a spectrum rich in lines until satisfactory 
results are obtained. This procedure takes a lot of time, especially if the rotation axis of 
the plate holder does not pass through the centre of the plate; moreover, it cannot guarantee 
that a further small change in one of the parameters would not yield a result still more satisfac- 


1 This most general case occurs e. g. for the ISP-67 spectrograph, which has four exchangeable cameras. 
Here, re-focusing has to be carried out not only for a new camera, but, even by changing spectral regions, as no 


single plate encompasses the entire visible spectrum. 
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tory. Thus it was felt that a method should be found for determining the optimum values 
of the variable parameters rapidly and with the least possible loss of time. 

The author's method is based on that of Hartmann's for investigating objectives (see, 
e.g. Martin 1950). Hartmann's method requires a diaphragm with two apertures to be placed 
behind the objective investigated. A sharp image of the distant point source is obtained 
only in the focal plane, whereas elsewhere two light spots result. Obviously, on placing 
such a diaphragm with two circular (Brode 1949, p. 301) or rectangular (James 1960, Arrak 
1959) apertures behind the spectrograph objective, sharp spectral lines will be obtained 
only in the focal plane; elsewhere, each line is splitted into two fringes whose mutual dis- 
tance b increases as defocusing e, i.e. the distance of the plane of observation from the focal 
plane, increases. By plotting the straight line b — (D), the value of lo, for which the spectrum 
is focused, provided / to be the only parameter, is obtained immediately. (This case is 
considered in papers quoted above). 

The disadvantage here consists in the fact that, as Í increases, the width of the fringes 
(which are wide enough by temselves) becomes greater, reducing the accuracy in determin- 
ing the quantity b. On the other hand, various spectroscopic investigations require a high 


degree of precision in focusing?. 


Principles of the method 


The present method consists in replacing the Hartmann type diaphragm by one of 
the kind shown in Fig. 4A (as decided upon after testing various other possibilities). This 
diaphragm contains two pairs of long slits (I) yielding two groups of equidistant Young 
fringes. The mutual distance of the groups increases with |e|. Between each pair of slits (I) 
there is a narrow horizontal slit (II) which serves for “marking” the central fringes of the 
group (see, Fig. 2). As the “marked” fringes differ from the remaining ones only by their 
intensity, they can be made sufficiently narrow by choosing the diaphragm appropriately 
(see the subsequent Sections). This circumstance, together with the fact that width of the 
fringes does not vary appreciably as [e| increases, provides the possibility of measuring the 
distance b between them (especially for large |e|) with an accuracy that is limited only by 
the grain size of the emulsion and by the efects of deformation of the latter during 
photographie processing? (see, e.g. Tolansky 1947, Ch. 14). 

The procedure to be adopted in measuring depends on the number of variable para- 
meters. 

A) If only one parameter e.g. l varies, one should make photographs containing a num- 
ber of spectra at different values of /,. Subsequent to measuring the distance between the 


* [f it is required to obtain e. 8. 90% of the resolving power of the spectrograph when using the so-called 
normal width of the slit as given by van Cittert criterion a, = Ak, the admissible defocusing eg, is given by 
legal = $ a„k (k denoting the reciprocal relative aperture of the collimator objective used) E. g. for a collimator of 
the spectrograph used in these experiments with k = 5, 54<eg<134 (according to the light wavelength), thus 
the requirements are very high indeed. With a wider slit the tolerance is greater but the resolving power lesser, 
as the normal slit yields maximum speed and resolving power. 

? At small |e| the groups overlap and, for certain values of e, the image is somewhat diffluent so that the 
accuracy of the reading of b can become somewhat lesser (see, e. g. Fig. 2b, photo for WW 
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Fig. 1. The straight lines b=b(/) for different values of p and different wavelengths. Six-slit diaphragm, 
3 ISP-67 spectrograph. 


“marked” fringes, b; with a comparator for each wavelength 4,, the straight lines b=b(/) 
should be plotted, which, if the other parameters are indeed fixed exactly, will then intersect 
n the /-axis in the single point required lọ. 

B) If two parameters are variable, e.g. l and g, several series of spectra should be 
aken at fixed values of y and variable / within each series. On obtaining the corresponding 
alues of b, the family of straight lines b=b(/) is plotted with different p as constant para- 

eters. These lines will intersect the /-axis at the points 4, ły, ls, ... (see, e.g. Fig. 1). By 
lotting g= (1) for different A,, a series of straight lines is obtained which, if the surface 
ontaining the foci for all A is plane or possesses the same curvature as the plate, will inter- 
sect at one point determining the lọ and gy required (see, Fig. 3)*. 
C) If the tree parameters /, L, g are variable, the procedure is again that of B), to be 
repeated. several times for different values of L. The results are plotted in the coordinate 
system 1, L, pọ, and the plane perpendicular to the L-axis is then determined whereon the 


4 In order to check the accuracy of the method, the values from five plates are shown in Fig. 3; usually, 
a single plate would have been sufficient. The presence of values for two wavelengths only is due to the fact that 
the remaining lines of the Hg spectrum appeared too weak on the orthochromatic plates, whereas it was consid- 
ered unnecessary to proceed to other measurements with a light source richer in lines, as the values of gy = 10.502* 
and /,—0.023 mm yielded excellent sharpness of the spectrum of the entire plate. These values were computed 
by least squares procedure for finding the equations of all the straight lines. The deviations from reproducibility 
in the values of pọ and I, thus obtained from the data of each plate separately amounted to + 0.02? and + 20m, 
respectively, which is five times better than the accuracy which can be determined by the spectrograph 


mechanisms. 
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Fig. 2. Example of interference patterns obtained with 6-slit diaphragm at different values of / and at constant p 
a) 4,—4358.34 A, i 
b) 4,—5460.73 A. Magnification about 26 times. 
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et] 
Fig. 3. Obtaining I, and q for 1SP-67 spectrograph. The points from five plates are marked e , O, x, D, V. 
straight lines p=g(/,) for the various wavelengths A, intersect in one point (the spectrum 
can then be focused sharper on the plane of the plate) or in several points situated as near 
as possible to one another (if the spectrum can be focused on the plane of the plate only 


approximately). This yields ł9, Lo and gy. 


Choice of the diaphragm 


When proceeding to determine the optimal values of the diaphragm parameters it 
should be held in mind that, as a matter of fact, the only really important fringes are the 
“marked” ones whereas all the remaining ones are a nuisance and therefore should be few 
and weak. Normally, at not too intense exposition, and utilizing a contrast developer, they 
cover but a part of the central maximum of diffraction pattern that would result from 
a single vertical slit I. Thus, the width of this maximum, as given by the relationship P— AF/yg 
(wherein F is the focal distance of the objective) should be small; this will be the case if 
the width 2y, of the slit I is large. Similarly, the distance between the slits I (3(72— 70) 
should be large, since then the distance between two interference fringes p —2AF/(n, — no) 
is small. However, with respect to the grain size of the plate, p should not be smaller than 
approximately 264. The two expressions jointly yield the number of fringes within the 
limits of the principal maximum, P/p=(n2—No)[2yo- Moreover, a high value of 2yo is con- 
venient with respect to the luminosity of the device (for this reason, the height of the slit I 
should exceed that of the prism); on the other hand, a high value of 2y, and i(gg —) 
would reduce adversely the depth of the area of interference given by the formula 
E150. 16AFP 2/yo(N2— No) (Emax denoted the highest admissible value of ¢). The foregoing 
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expression was obtained on considering the intersection of two diagrams of isophotes of 
diffraction patterns outside the focus, as resulting from either slit I separately. The numer- 
ical coefficient 0.16 is an empirical one.? 

The horizontal slit II should obviously possess a horizontal edge of length 7j, — nọ in 
order that the light emerging therefrom may modulate, say, three interference fringes. Its 
vertical dimension, which determines the degree of modulation should be chosen in such 
a manner that the area of the slit II shall constitute 1095 of that of slit I (as is found by 
trials). 

The distance between either group of slits should be large in order to obtain high sen- 
sitivity of the device (the variation of b corresponding to a given change in e increases with 
the distance separating the slits). The diaphragm, which is best made of metal, ought to be 
set up in the spectrograph so that the optical axis of the latter shall pass through its centre 
and that the interference fringes from either group of slits observed under strong magnifica- 
tion shall be as sharp as possible. The sharpness is regulated by rotating each group of 
slits about the horizontal axis. It is especially important to have the possibility of such ro- 
tation if the spectrograph slit is curved (as was the case of the ZSP-67 spectrograph). 

The light source should lie exactly on the axis of the spectrograph, as even a slight 
deviation causes considerable differences of the light intensity of either group of fringes. 
This effect (already mentioned by James (1960)) can serve for setting up light sources on 
the axis to a high degree of precision; this, in turn, is known to affect strongly the sharpness 
of the lines. 

By using a light source rich in lines, the entire focal plane can be found with greater, 
aecuracy. The lines, however, should not be too near one another, as then their interference 
patterns might overlap, rendering measurement more difficult. 

As the distance between the fringes p depends on the wavelength, the method is best 
adapted to the case when the plate encompasses a not too great region of the spectrum, 
(e.g. the visible spectrum or part of it). If this region is large (e.g. from 2000 Á to 7000 Á), 
a single diaphragm will not yield sufficiently dense fringes throughout the entire region 
unless an especially fine grain plate is used. It is then necessary to utilize successively two 
different diaphragms. The suitable diaphragms can be more easily chosen for spectrographs 
with large prisms. 

In addition to focusing spectrographs, the present method can be used for finding the 
focal lengths of lenses (indeed, this was done with entirely satisfactory results); however, 
here, the accuracy cannot compete with that of the well-known three-slits method of Zernike 
(1950) which has been used for this purpose e.g. by Renét and Vićnot (1954). 


APPENDIX 


Computation of the intensity distribution throughout the interference pattern, 
In order to evaluate the intensity distribution at an arbitrary point P of the focal plane 77 
(Fig. 4A) in the case of a point source S and an ideal objective (of focal length R equal to 


5 If the light is not perpendicularly 


incident on the plate, the sine of the angle of incidence should appear 
in the formulae. 
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the segment OC) covered by a diaphragm with six slits distributed as above in the plane 
7'B', a method was used resembling the one proposed by Françon (1956, p. 292) and applied 
to Fraunhofer diffraction as occurring in the present case. 

/ For the amplitude of the oscillations at point P, the following expression is obtained: 
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wherein K—2z/A, and all other notations$ are those of Fig 4A. On integration and 
squaring, Eq. (1) yields the intensity at P: 
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In order to evaluate the intensity distribution at points P’ at the distance e from the focal 
plane IT, it is required to add” to the exponent of each integral in (1) dependent on f the 
quantity 4Keß’?, and to the exponent of those dependent on y’ — the quantity $Key 

In order to integrate the new expression, we transform all integrals to the Fresnel type 


F(w) = fow (: 3 e) dr by introducing the new variables t}=f'"—y'/e and Ty=y —z'[e 


0 
and using substitutions similar to those used by Sommerfeld (1959, p. 217). Finally, Ep. 


is obtained in the form 


Ep. = ree) — F(w))] Y [FEw) — F 1j] + 


Š 2 
+ [F Cw) — FQwg) > [F — Few] (3) 


6 All primed quantites are in angular measure, e. g. 2y5 is the UE. width of the slits I. The same quan- 
tites without primes in the ensuing formulas are expressed in linear measure, thus e. g. yo = Yol m Po = Bln etc. 


The quantites z”, y”, €, are in linear measure. 
7 This substitution results from Eq. (67.6) (Françon 1956, p. 271) car z = p°? + y^. 
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wherein: 
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Eqs. (2) and (3) are considerably simplified on putting y' —0, i.e. it we consider the intensity 
AA e sin f(y’) 
distribution along the z'-axis. Indeed, in Eq. (2) two expressions of ihe type — —-—— now 


fO 


vanish, whereas in Eq. (3) each difference of two Fresnel integrals preceding the sign >; 
reduces to a single integral, as then ?w4— —!w,, ?wg— —!wg and F(w) = —F(—w). 

In Fig. 4B, the intensity distributions in the three planes e —0, e —20004 and & —5000A 
are given as evaluated by substitution of the given numerical values in (2) and (3), at y’ =0.8 

For the case e=0 (2), the principal minima are obtained when Az'yy/R=k,r and 
simultaneously Kz'ygya| R =k, wherein ki, I = +1, +2, +3, ... (if a is an integer, this 
occurs beginning by k =1). The secondary minima occur where one cosine in the first 
term of (2) and the cosine of the second term vanish simultaneously. For the case con- 
sidered this occurs for: 


ER = 100004,/14 J=+1, +2, +3... (principal minima) 
z, = 1002k/14 k= +], +3, +5... (secondary minima) 


For 40 (3), the positions of the minima cannot be predicted but had to be found 
by numerical computation, as also the other points of the curves. From Fig. 4B it is seen 
that, notwithstanding the totally different formulae, the intensity distribution exhibits the 
same trend for £=0 and for ¢40 and the minima occur at the same intervals. 

The perceptible secondary maxima can be easily grouped by five (the central maximum, 
the two neighbouring weaker ones and the two quite weak outer ones). These groups are 
separated by a distance of 10002/14 which is the distance p between the Young fringes that 
would result if, instead of six slits there would be only a pair of slits I (with the same distance 
in between). Thus, each fringe on the plate should be expected to present the above struc- 
ture, which, however, fails to appear, as the optical system is no ideal one, because of the 
finite dimensions of the slit and the too small resolving power of the emulsion. In the region 
considered, as seen from Fig. 4, the strongest groups of maxima correspond to the “marked” 


* These values (not greatly differing from the real ones) make it possible to simplify further Eq. (2), as 
yielding 2%, = 1 + Nz, thus allowing one cosine to be placed before the brackets. 
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fringes. As |e| increases, so do the distances between them, shown in Fig. 4 as b, and bs. 
pvo. £M cr : f tte; 
For higher |e|, as estimated but not shown on the graph, these fringes become weaker than 
i inni i 7 s i sitions can be just de- 
their neighbours beginning from a given value of e, so that their positions ca j 
termined by those of these neighbours (this can be seen e.g. in Fig. 2). a 
š ; USE os ? ; RSI -—— 
Stricly speaking, we are dealing in experiments with the illuminated slit, and o 
a point source of light. Thus, in order to obtain expressions corresponding to this case, 


it js necessary to evaluate the integrals 


, 


50 Yo 
[Erd and f Epdy 


-» -5 
(Françon 1956, p. 284), where yọ is half the angular height of the slit. The integration of 


the first integral yields 


2 
Ep = UE (& B?b + éC? + 25 Bc) (4) 


with B denoting the first and C the second term of the expression in the brackets of formula 
(2) if y' «0 is assumed therein, and 
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For sufficiently large? yọ, compared with A, we have 5, =£,=ë,=1. We see that only if 
b->1, will (4) be equal to (2), to within a constant factor before the brackets. However, 
since in fact b differs from 1 quite considerably (in our case b=0.02), the trend of the cor- 
responding graph (Fig. 4B, the curve denoted s=0, yg=1 mm) will undergo modification. 
Indeed, the “marked?” maximum is seen to be increased considerably as compared to the 
remaining ones (in agreement with observation). Other features, such as the shift in the centre 
of gravity of fringes 2 and 3 by 1002/14, or the decrease in intensity of fringe 3 with respect 
to fringe 4 could not manifest themselves on the plates, as the very structure of the fringes 
is unresolved (see above). 

It would seem that integration of (3) should yield similar modifications in the remain- 
ing curves (as inspection of Phot. 2 may suggest). However, we refrained from computa- 
tions concerned as this would involve integrals of a product of Fresnel integrals, and even 
with approximate formulae no simple expressions would be obtained (these would contain 
a linear combination of trigonometric function and Si and Ci). 


° For the diaphragm used, already at yy = 1 mm, ë, = 0.999493, £, = 0.974749, £, = 0.999994. 
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Maybe it is worth mentioning that, with the slit II not at a middle distance between 
the pair of slits I, but e.g. placed asymmetrically nearer to the outward slit I, the interference 
pattern would undergo deformation (the outer fringes of a group would then present an 


. intensity comparable to the central ones, and those contiguous on the latter would become 
quite weak). 
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On the assumption of small rclative variations of the vessel's radius, the amount of the 
fluid passing in unit time is found to be given by the telegraph, equation. Moreover, this equa- 
tion is found to be fulfilled, in the foregoing approximation, by the radius of the vessel as a func- 


tion of the time and space coordinates. 


I, Introduction 


The problem of non-stationary flow of a viscous fluid through vessels possessing elastic 
walls has been dealt with by a number of authors from the point of view of physics. The 
problem, too, is of high practical importance for the dynamics of arterial circulation. The 
equation of flow for the case of an ideal fluid and on the assumption of very small relative 
variations of the radius of the vessel is identical with the wave equation. It was derived by 
Korteweg (1878) and King (1947). It leads to the formula for the propagation velocity of 
a disturbance. In the case of blood circulation in the arterial system, this velocity is that 


, where c 


of the pulse wave. The velocity is given by the following formula: ¢ = Vi 5 
ELS 


is the velocity of propagation of the disturbance, E —Young's modulus for the wall of the vessel, 
R, — the radius of the vessel, ó — the thickness of the wall, and o — the density of the 
fluid. The starting point for our consideration of viscous flow was provided by the Navier- 
Stokes equations. The solution of the problem of flow through vessels having rigid walls, 
on the assumption that the pressure gradient and fluid velocity are Ae” type functions 
of the time, is due to Lambossy (1952) and Womersley (1955). For the case of elastic vessels, 
the same problem has been considered by Morgan and Kiely (1954) and by Womersley. 
The above-mentioned papers deal with linear theories. The non-linear case for an ideal 
fluid was dealt with by Lambert (1958) and, for a viscous fluid, by Miękisz (1961). The 
equation of the problem under consideration is derived in the last-mentioned paper quite 
(821) 
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generally from a variational principle. The equation is of the form 
gener: 


TR) 320 Q? | 
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where Q is the mass of the fluid which crossed the section x=const. of the vessel while the 
radius increased from Ry to R(x, t); Rọ denotes the smallest value of the vessel’s radius, 
b — a quantity related to the velocity profile (b—1, for a parabolic profile b —4/3 and for 
an ideal fluid 6=1), whereas 7 is the viscosity coefficient of the fluid. The quantities E, ó 
and o are the same as in the expression for the velocity of propagation of a disturbance. 


II. Hydrodynamic-electric analogies 


Eq. (1) is a quasi-linear partial differential equation of the second order, of the hyper- 
bolic type. The velocity of propagation of a disturbance is now given by the equations of 
the characteristics, which are of the form 


=) Ó : 
= ance z (2) 
[S 12 meer ic Sy 

age = 0, 


It is easily verified that, on the assumption of Q. < Rizo this yields 


dx dw BEAT 
(ela vom Kaci ©) 


-is the average velocity of the fluid (dS is a surface element, and S — thearea 


: vU in 
where v J 


of the sm: Womersley proved that the ratio v/c is of the same order as the ratio 


r/R, where R denotes the radius of the vessel, and r — its variation. The assumption of 
0O < R? no is one le ading to the linear approximation, where the equation of our problem 
CH 0 70 PI p 


isa linear one of the type 


30, 8y AQ Eó 92*Q 0 4 
02 BOORI Ot — 2boRQ,O0x* — (4) 


This yields the propagation velocity of a disturbance in the form c — _Bó_ ‚which trans- 
2boR, 

forms to Korteweg’s ideal fluid formula at b=1. 

On differentiating Eq. (4) with respect to the time, an identical equation is obtained 


for the amount of the fluid passing in unit time, which we shall denote by W. Thus 


WV Mm OW  Eó OM enn i 
9 ^ boRg 9t 2beR, 952 — 6) 
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Eq. (5) is the telegraph equation, known from the theory of non-stationary electric currents 
(Simonyi 1956). Hence, the well-known analogies between motion of a viscous fluid and 
flow of electric current can be generalized for our problem of viscous flow through elastic- 
walled vessels. In the theory of non-stationary electric currents, the current as function of 
the coordinates of the position and the time fulfills the following equation: 


9" Roi 1 00% 
oo Iota Lil 2%ż 


=0 (6) 


with Z, R and C denoting, respectively, the densities of the self-induction, resistance and 
capacity. Analogously introducing the notions of self-induction, resistance and capacity, 
L*, R*, C* in the present problem, and on equating the coefficients of Eq. (6) and Eq. (5), 
we have 


2 
co id, eC że 


Ir ~ boR3 EO (9 


The idea of generalizing the hydrodynamic-electric analogy to viscous flow through elastic- 
| walled vessels, which led to the telegraph equation, was put forward by R. Ingarden in 
1957 in a discussion with the present author. Independently, these analogies have been 


also considered by Ródenbeck (1960), Beier (1960) and Taylor (1957). 


IIT. Dependence of the radius of the vessel on the coordinates x, t 


The problem of how the vessel’s radius depends on the coordinate of position in the 
case of stationary viscous flow was discussed by Rashevsky (1945). The method utilized in 
deriving the equation of non-stationary flow also leads to the solution of the problem of ' 
the dependence of the radius on both coordinates in the general case considered. 

The mass of the fluid present within an element of the vessel of length dx is 


90 


ax 


R=R, V 1— gum (9) 


On expanding the right-hand side of (9) and on assuming Q,— Reno, we have, in the first 


approximation, 


R?zodx Ri mods — dx (8) 


whence 


= Qx 10 
a 2Ro uo 
By (10), differentiation of Eq. (4) with respect to x gives, for the radius of the vessel, 


the following equation: 
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Introducing here for convenience the quantites L*, R*, C*, we have 


OER pid? aR WERL (12) 
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We shall restrict ourselves to solutions of Eq. (12) of the type 
R = A (xe (13) 


where 4(x) is independent of the time. This is equivalent to assuming that the function 
representing the dependence of the radius on the time is a periodic function at its boundary 
and can be expanded in a Fourier series. This is indeed the case i.a. in applying our method 
to blood circulation, where the rhythmic pumping of blood into the aorta by the heart 
gives rise to a periodic dependence of the variation of the vessels radius at its boundary. 


Substitution of Eq. (13) in (12) yields 


TA — y AG) (14) 
with 
y? = — o?*I*C* + io R*C* 
whence 
y = + Vio R*C* — o*L*C* (15) 
The general solution of Eq. (14) can be written as follows: 
A(x) = R, e^"* + Ri e™ (16) 
Putting y = + (z + if), we have 
A(x) = Re 2-188 4. RV eax Hia (17) 


The quantity & is the damping constant, and £ is related to the phase velocity of propagation 
of the disturbance as follows: 
sz api (18) 


Introducing (17) in Eq. (13), and taking account of (18), we have the general solution, which 
can be written thus 


x 


R-—R,e-**e^ ( >) + Ri eos e” ( Ą >) (19) 


By Eqs. (16) and (7), and on simple computations, the damping constant « is obtained in 
the first approximation in the form 


with » — 21 and GEE Łó = 
0 2boR, 


* = Ribc (28 


o?w0?b? Š 


Eq. (20) holds only for vessels fulfilline the condition Ry = Ver 
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In the case of 


motion of blood in the vessels, where 7 =0.04 g/cm sec, 9 =1.02 g/cm? and o =2z 1.16 Bec, 


we have the following condition for Ry: Ry>0.18 cm. The approximate formula (20) can 


be applied in blood dynamics in the case of the large vessels only. For the narrow vessels, 


the coefficient should be computed from the exact formula 


Qe y -1otrcu 


T 
2 


V (R?*2 + 92152) 2*2 
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ON THE DISPERSION OF THE PROPAGATION VELOCITY OF 
DISTURBANCE IN A LIQUID FLOWING THROUGH ELASTIC. 
WALLED VESSELS 


By S. MIĘKIsz 
Chair of Physics, High School of Agriculture, Wrocław 
"(Received April 24, 1961) 


From the equations of non-stationary viscous flow through vessels having elastic walls, 
| a formula for the velocity of propagation of a disturbance is derived. The velocity is found to 
depend on the frequency. Dispersion of the velocity is due to the viscosity of the fluid. A formula 

for the group velocity is also derived. 


I. Introduction 
The equation of non-stationary viscous flow through vessels possessing elastic walls is, 
in the linear approach, of the form (Miękisz, 1961) 
aL vi nao] "BOO". 1 
0 * beRą 0: 2boR, dx? (1) 


where JV is the mass of the fluid passing a given point in the vessel in unit time, 7 — the 
viscosity coefficient of the fluid, o — its density, Æ — Young's modulus for the walls of the 
vessel, ó — the thickness of the wall, R — the radius of the vessel and b — a constant 
dependent on the velocity profile. The foregoing is identical with the wave equation and 
is known as the telegraph equation in the theory of non-stationary electric currents (Si- 
monyi 1956). This provided the possibility, in an earlier paper, of proposing hydrodynamic- 
slectric analogies for the non-stationary flow of a viscous liquid through vessels having 
elastic walls, thus introducing the notions of self-induction, resistance and capacity in the 
problem under consideration. By comparing the coefficients in the equation of either pro- 


sess, the following formulae result: 


Rx 85 


2boR 
ct rU TR (2) 
Dx bok? Eö 
vhere R*, L* and C* are the linear densities of the resistance, self-induction and capacity, 
espectively. Hence, some relationships of interest, appearing in the theory of non-stationary 
urrents, can be transferred to the present problem. 


(827) 
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II. Dispersion of the velocity of propagation of a disturbance 
On writing the solution of Eq. (1) as follows: 
= Jue "R 


y is readily seen to be given by the formula 


y= + |; TO LRG = TORES pe + (a + ip) 


By simple computations, ff is obtained in the form 


C ——— TEŻ ; 
p= V3 o? L*C* + = y R*? + wils (4) 


The quantity f is generally known to be related to the phase velocity of propagation of 


a disturbance by the relationship 


j 1 
(5) 
O l 
Eq. (3) for y is easily rewritten thus: 
w VIF aiii 
y = io J LG ( sp) (6) 


We shall restrict further consideration to the case of weak damping, when R*/@Ll*<1. 
Expanding the right hand side of (6) in a power series, B is approximated by 


Reed 1 Ra 
= / IxC* zr 
B= o / LX ( = za) (7) 
By Eq. (5), the phase velocity is obtained in the form 
l UERS 
A > di] MAP 
i x 80? r | (8) 


Introducing the expressions of Eq. (2) into Eq. (8), we finally obtain the following expression 
for the phase velocity of propagation of the disturbance: 


8»? 
Vibe ( 7 iis] (9) 


Pó 
wherein c "Va RC As the phase velocity v is seen to be a function of the angular 
400g 


velocity o, we are dealing with dispersion of the velocity. It is the viscosity of the liquid 
that accounts for this effect. The group velocity is known to be given by 


_ do 
^. dk 


u 
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Moreover, we have 


4 
u == 
1-0 de = 
v do 
In the present approximation, substitution of (9) in Eq. (10) yields 
A, 8»? 
u= EMM Ae Rabo? (11) 


III. Group velocity and phase velocity of a pulse wave 


| The investigation of non-stationary flow of a viscous liquid through vessels possessing 
elastic walls presents considerable practical interest with respect to the dynamics of blood 
circulation. Here, the velocity of propagation of the disturbance is that of the pulse wave. 
The proof of the existance of velocity dispersion was based on the assumption that R*/@L* < 1 
Now, by substituting for R*/L* the value of this ratio from Eq. (2), we can obtain the con- 
dition to be fulfilled by the radius of the vessel. This condition is of the form: 


85 
R ç 
o> V; 2 (12) 


On substituting in (12) the numerical values characterizing the dynamics of blood: 
n=0.04 g/cm sec, 9=1.02 g/cm?, b=4/3 and w=2an, where n=7/6 is the fundamental 
frequency, we obtain R,<0.18 cm. Thus, the foregoing consideration concern large vessels 


ME š : : w ; 
only. On substituting the dimensionless quantity x = R] / — , as introduced by Womer- 
y 


sley(1955) in haemodynamical theory, into Eq. (9) determining the phase velocity, we have 


8 
je E = =) (13) 


The numerical values of the ratio v/c for different values of œ are assembled in the following 


Table: 


| | 
a | 1.8 | 2.0 2.2 | 2.4 | 2.6 | 28 | 3.0 3.9 | 4.0 | 5.0 
f | ` 
ać | 0.57 0.71 0.81 0.86 0.9 0.93 | 0.95 | 0.97 0.98 0.99 
Ś | d | 
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Fig. 1 shows the dependence of the ratio v/c on a according to the equation (13). Fig.2 i 
that of the dependence of the ratio v/cg on z, according to the results of Womersley (1957 


Eô Bo u 
(a = V six; and c — V si = 0.86 co) ° 
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ON THE ELECTROMAGNETIC POTENTIALS OF AN ELECTRO- 
MAGNETIC DIPOLE IN MOTION 


By A. BIAŁAS 


Institute of Physies, Jagellonian University, Cracow 


(Received April 26, 1961) 


Starting from the Lorentz solutions of Maxwell's equations, expressions are derived for the 
electromagnetic potentials of an electromagnetic dipole in arbitrary motion. By utilizing Dirac's 
0—function in the description of the field sources, these expressions are generalized to multipoles 
of arbitrary order. For a multipole at rest: they become identical in the wave zone with the well- 
known formulas for multipole radiation despite the fact that the time dependence of the multipole 


moments is widely arbitrary. 


l. Introduction 


The well-known Liénard-Wiechert formulas for the electromagnetic potentials of 


a charge in motion 


ee ee (1.1) 


fail to take into account its structure. Hence they can be considered to represent the zero 
approximation neglecting all terms dependent on the dimensions and internal motion of 
the body under consideration. That is why it was felt to be of interest to derive formulas 
analogical to (1.1) and taking into account the structure of the field sources. The present 
paper is devoted to this problem. Equations for the electromagnetic potentials of a system 
of charges and currents in motion will be derived in the first approximation, t.e. taking 
into account the terms linear in the dimensions of the system. A body whose description thus 
given is sufficiently strict and whose total charge vanishes will be termed an electromagnetic 
dipole. The method to be used here is a generalisation of the well-known method used in 
deriving the formalas of (1.1) (see, e.g. Lorentz 1909). The same expressions are subsequently 
derived utilizing Dirac’s ö-function for describing the sources of the field, making it possible 
to generalize our formulas to multipoles of an arbitrary order. The last Section is devoted 
to the relationship between the formulas obtained and multipole radiation. 
(831) 
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We now proceed tc a somewhat more strict formulation of the assumptions: | 

1. The ratio b/R, wherein b is the diameter of the body, is considered to be an infinitely 
small quantity of the first order; the calculations are carried out to the first approximation. 

2. The following quantities will also be considered to be at the most infinitely small 
of order one: 

a) the ratio of the relative velocity of any two points of the body and the velocity of 
hight; 

b) vb/c, where v' is the acceleration of any point of the body. 

If a system of charges and currents moves with constant velocity along a straight line, 
then the current density can be relativistically-invariantly represented as the sum of a part 
"transferred" and of a part conducted” (see, e.g. Becker 1936). Denoting the total current 
density by j and that of the current conducted by Jj. we have 


jJ jr exw | (1.2) 
wherein og is the real charge density, and 
w — the velocity of motion of the body. Similarly, we have 
ire w 
et (1.3) 
a 


with o denoting the total charge density. 

If motion is not rectilinear, Eq. (1.3) ceases to be valid. This becomes clear at once 
if we consider that then conservation of charge is not generally applicable. Hence, the 
expression of Eq. (1.3) is to be taken as an approximation of the corresponding strict ex- 
pression. Although we have no knowledge of the exact form of the latter, it suffices for 
the present purposes that any correction to be added to (1.3). by the foregoing assumptions, 
is infinitely small of the first order. Thus, we have 


CPE Or (1.4) 


l 
QL = -y JL: t (1 B) (1.5) 


with € — the velocity of any point of the body, and f — an infinitely small quantity of order 
one. The construction of (1.5) should be such as to ensure the fulfillment of the continuity 
equation by jj and oj: 

JoL 


o * dvjL=0 (1.6) 


(obviously, we require no more than that (1.6) should be fulfilled with first order accuracy). 
By (1.6), we have 


Jed =0 (1.6a) 


With respect to the separation of (1.2) and (1.3), and taking into account assumption. 2a, 
we can also split any electromagnetic dipole into an electric and a magnetic part. If jr=Q 
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we say that we are dealing with an electric dipole. If 94, —0, we term it a magnetic dipole. 
Thus, the condition for an electrie dipole is | 


j=w (1.7) 


and, for a magnetic dipole, 


o =~ (L+ f) (1.8) 


2. Derivation of the Liénard-Wiechert potentials by Lorentz’s method 


In order to make the reader acquainted with the method to be followed, we now proceed 
to derive Eqs. (1.1) according to Lorentz (Lorentz 1909). We begin by the well-known ex- 
pressions yielding the solutions of Maxwell's equations through retarded potentials! : 


g(XYZT) = | ELLOS dV’ (2.1) 
acerzn | ERA zy. (2.2) 


wherein the following notation has been adopted: 
RIR R(X- Voy, X 


X Y ZT — the point-moment at which the field is being determined, 
xy z  — the current point of the integral. 
The moment 7 is determined from the equation 


CTR (2.3) 


and will be henceforth termed the “effective moment”. 
As now (1.7) holds, Eq. (2.2) can be rewritten thus 


WAŻ d wat) dV’ (2.4) 


A(XYZT) = ale 


In order to obtain the expressions required, it is necessary to carry out the integration of 
the right hand sides of Eqs. (2.1) and (2.4). As, in general, to each curreni point of the 
integral corresponds a distinct active moment, this cannot be effected directly. We have 
first to choose new variables in such a manner that all subintegral expressions shall correspond 
to one and the same moment, independently of the point of integration. 

We arbitrarily determine a point M within the body (henceforth to be referred to as the 
centre). Let us assume (this does affect generality) that, at the moment 9, the centre is al 
the origin of the coordinate system and that 0 is its effective moment. This yields 


R=T-9 (2.5) 


1 Henceforth we put c— 1. 


834 


As M is the origin, we have 


R=YX?+Y2+Z2 R=(X,Y,Z) (2.6) 


Consider now an arbitrary neighbouring point M’ of the body. Let us denote its effective 


moment by 7. The equation for 7 is of the form | 
IR —r()|-T—« (2.7) 


Here, r(t) is the radius vector (position vector) of M’ at the moment r. Thus, obviously, 


the transformation 


r(r)j=r+vdt | At —7—0 (2.8) 


yields the first approximation of the effective position of the point M' as given by its position 
r at the moment 02. As it is our intention to neglect the dimensions of the body in our 
ultimate result, we can indeed restrict ourselves to the first approximation in (2.8). This 
is so because only terms of the first approximation can contribute finitely to the transforma- 
tion Jacobian which we have to set up in order to transform the volume element dV’ = 
=dx (v)dy' (r)dz' (v) to dV —dx dy dz. In this approximation, simple computations yield 


R-r R-v 
IR — r(?)| = R — Xd ds : At (2.9) 
By Eqs. (2.5), (2.7) and (2.8b), we have now 
R-r 
= 2 
It hic hts (2.10) 
We now introduce the notation 
R—R.v=!] (2.11) 


With (2.10) and (2.11), computation of the Jacobian of the transformation (2.8) yields 


J = det E) aie (2.12) 


0 XR 


On substituting this in the integrals of (2.1) and (2.4), and taking account of 
Joy) dV = e (2.13) 


wherein e is the charge of the body, we have with an accuracy to terms of order zero. 


NM ev e 
i) beets ił (2.14) 


hese are precisely the expressions we required. 


Henceforth if a quantity is not provided with an index showing to which moment it is related it will be 
meant to be related to the moment 0. 
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The foregoing formulas can be rewritten in a way to lay stress on their invariant charac- 
ter?: 
eu” 
ou 
wherein u” denotes the four-velocity of the body, and (X, Y, Z, T— 0) =R” — the four-vector 
connecting the effective point-moment of the body and the point-moment at which the 


A^ (2.15) 


potential is to be evaluated. By (2.5), the four-vector R* fulfills the equation 
IN, JC FO (2.16) 
As we consider only retarded potentials, the following condition holds too: 


R,>O (2.17) 


3. Change of variables in the first approximation 


In order to determine the electromagnetic potentials of the dipole, we once more make 
use of Eqs. (2.1) and (2.2). The considerations will be fully analogical to those of the preced- 
ing Section, the sole difference being that we now proceed in the first approximation. Thus, 
we have primarily to find a transformation reducing the subintegral expressions to a single 
moment. This transformation must be found in the second approximation, so that we can 
compute its Jacobian in the first approximation. With the notation already used, we have 


rf) =r + VAT +r'AT + Š (AT)? (3.1) 


Herein, ** denotes the velocity of the point M” with respect to the centre of the dipole at the , 
moment 9, and a is the acceleration of the center of the dipole. The time interval AT is 


given by y 
AT—7—0 (3.2) 
AT differs from At of the preceding Section by an infinitely small quantity of order two: 
AREA (3.3) 


On substitution of (3.3) in. Eq. (3.1), and taking into account that, r is a small quantity 


of order one, we have, to the second order, 
a 
r(r) = r+ At + pa + At + 3 (At)? (3.4) 


Employing the foregoing formulas, we compute & quite similarly as in Section 2, and, 
subsequently, the expression appearing in the integrals of (2.1) and (2.2): 

E ur is Dart m Hr tier) (3.5) 
|R — r(z)| pe is la i? Í 

Hitherto, the calculations could be carried out for the electric and magnetic dipoles jointly. 
We shall now proceed to deal with either of these two cases separately. 


3 Here and henceforth we employ the metric — — — +. 
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1. Electromagnetic potentials of an electric dipole 


In the case of an electric dipole, as shown in Section 1, the current density is related 
to the charge density by Eq. (1.7). Hence, the potentials can be determined by means of 
Eqs. (2.1) and (2.4). To this purpose, we carry out the transformation (2.8) in the functions e 


and @ appearing in those formulas. We have, as is well known 


D 
o(r + vAt, 0 + At) = o(r, 0) + T At (4.1) 
Here, Dr denotes the substantial derivative. As o fulfills the equation of continuity 
D 
23 + odivr"=0 (4.2) 


we have 


o div (R - r) 
] l 


o(r + vAt, 0+ At) = o(r, 0) (4.3) 


Substituting this in (2.1) and invoking (3.5), we obtain 


1 221—324 (R - kt ad ja : 
C= E 0)d V + : E a) I : r)od V — p L : ve)od V 4 A (R:r)odl_ (4.4) 


We now introduce the notation 
d = f ro(r, Yav | (A.5) 
In what follows, we shall term d the electric moment of the dipole. With respect to Eq. 
(4.2), 
d' = | r'efr. ar (4.6) 


where denotes the differentiation with respect to the time. From the foregoing, and re- 
membering that e=0, we have 


l-v+R-a v:d R-d 
p= — p — Rd — + (4.7) 
We now evaluate 4. By 
j à Do 
v(r + vdt, 0 + At) = v(0, 0) + 1'(0) + Dr At, (4.8) 


we obtain. in full analogy to the earlier result, 


A=vo+a Br + = (4.9) 


Eqs. (4.7) and (4.9) provide the required expressions for the electromagnetie potentials of 
an electric dipole. 
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5. Electromagnetic potentials of a magnetic dipole 


A magnetic dipole is characterized by the fact that the current density within it is related 
‚to the charge density by Eq. (1.8). Thus, Eq. (2. 1) now assumes the form 


Qi 2 T 0 s. er A 
A E Y, 2,2) v G^ y zv) (+ a) y, 


R qa (5.1) 


which, together with (2.2), will serve for determining the potentials required. To this pur- 
pose, we change the variables in the function j(«’, y'; z', z): 


j(r + v At, 0 + At) = j (r, 0) + E: At (5.2) ` 


Substitution of (3.5), (4.8) and (5.2) in Eqs. (5.1) and (2.2) yields 


INCUN RR: 
Ant | iar =a n SL ZYD (5.3) 


IE UR: 
> u = fe: GR - r) = £ ZI DR yar (5.4) 


In (5.4), the relationship (1.6 a) has already been taken into account. 


We now proceed to transform the expressions of (5.3) and (5.4). For this, we shall 


compute 
ip (r: RyjdV/ (5.5) 
Note that, in the system at rest i.e. such that p=0, we have (Średniawa 1960) 
f iar — —4Rx f (rx fav (5.6) 


The vector 4 f r — jdV as computed in the system at rest will be denoted by w. We shall 


refer to it as the magnetic moment of the dipole in the resting system. All quantities referr- 
ing to the rest system will henceforth be given the lower index ç. Thus, 


f R: rb dV, = u xR (5.7) 


In order to evaluate the integral (5.5) in the system in which the dipole moves with veloc- 
ity e, we express the quantities r, j, dV by means of ro, jo, dro. In accordance with the 


transformation rules, we have 


j=lt 6 00-0 (5.8) 
I š: 

a=To— „5 (ro^) jA (5.9) 

ay = 0 | (5.10) 


Ya 
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with 
l 
I ox kil 
y y1 — y? ( ) 
Additionally, we introduce the vector m related to w as follows: 
45 Ko, 1 N 
m = uw — = (e * 1) ene (5.12) 


Substitution of (5.8), (5.9) and (5.10) in Eq. (5.5), together with Eq. (5.12). readily yields 
the following relatioship: 


[ir aV — mx R+ (R- v)(e x m) (5.13) 
From (5.13), it follows immediately that 
[we DR nar = v f jR- r) dV = v: (m x B) (5.14) 
Let us also compute if jdV. For this, we make use of 
[dav = f v div jay (5.15) 


The foregoing relationship holds for an arbitrary vector j vanishing outside the region of 
integration. With respect to the equation of continuity, we have. to the first order, 


a) ` 

fr div j dV = [ro (j-v)dV = - ra - p) dV (5.16) 

Considerations quite similar to those leading to the evaluation of the integral (5.5) yield 
f r(j.v)dV -vx m (5.17) 

Thus, 
MEG. 
JEJ > (v x m) = (vx m) (5.18) 
( 


By means of Eqs. (5.13), (5.14) and (5.17), we are now able to rewrite (5.3) and (5.4) as 


follows: 


Wu m T m 
A = — basi x R+ R(vx m)] + D à R PA. 


a4 9 PO R: (m x v) 
Dil 


(5.19) 


6. Alternative method of deriving the formulas for the electromagnetic potentials of 
a dipole in motion 


It is often convenient to deal with an electromagnetic dipole as if it were a singularity 
of the sources of the field. In this case it is posible to use for its description Dirac’s ó function. 
Obviously, this method cannot be said to be quite exact, as we have no full knowledge of 
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what transition to the limit we are effecting here. Notwithstanding, it is in general use e.g. 

as applied to the gravitational field (Infeld 1954, Tulczyjew 1959), and is felt to present 

the advantage of simplicity. In addition, it yields the result in relativistically invariant form. 
Green's function for the Maxwell equations | 


LI4”=4aj* (D=9,9% (6.1) 
and corresponding to retarded potentials is given by (see, e.g. Thirring 1958) 
l Sag z 
Dret(x#) = 157 (“PX u) for x4 >0 (6.2) 
0 for x4 <0 


Thus, the solution of Eq. (6.1) is of the form 


A" (x^) = 4n li D SE IgA) da (6.3) 


It is our aim to find these solutions in explicite form. The field source is an electromagnetic 
dipole. We shall assume that the latter is described by a four-vector of current of the form 


j“ = | dsD"( 9,8" — 2"(s)) (6.4) 


wherein z^ —z"(s) is the equation of the world line of the dipole, and D"'(s) is an antisym- 
metric tensor. D"' will be termed the tensor of moments. 

The choice of the form (6.4) for the current comprises both the electric and magnetic 
dipoles. The kind of dipole depends on the form of the tensor D“”. In the case of an electric 
dipole, D"" in the system at rest is of the form 

MEET 
0 0 Wór 


[DRO == 6.5 
ç ny ape ee) 


Sse peo 


Hence, in an arbitrary coordinate system, it can be represented thus 


Dia Du = Duż (6.6) 
where D" is a four-vector whose components in the system at rest are: 
D$ — (Dy, Dy., Dy”, 0) (6.7) 
D^ is seen to fulfill the equality 
D'u,-— 0 (6.8) 


For a magnetic dipole, D'" in the rest system is of the form 
W MO Ms 0 
ue | eo MBE 0 re 
eye sooo | 


0 Du 00 
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This condition can be put in invariant form: 
Mu, =0 (6.10) 
Having defined the four-vector of current density, we can proceed to solve Eqs. (6.1). 
wilka : BT NE © ; 
This consists in carrying out the integration in (6.3) after replacing therein j“ by (6.4). By 
the well-known property of the derivative of the 6-function: 
| f(x) 969 dx = — f' (0). (6.11) 
we have 


A mag JD") 2, pre (xf 2"(s)) ds (6.12) 


On differentiation and integration per partes, this yields 


ES | 
ŻE ne [Dret (x^ — z^ (s))] E (x^ — z^ e» | 2(x, — z, (s)) D" | = 

d [De — z,(5)) | 6.14 

A ret 3 
B IE ds | (x, — z, (s)) w _ (AA 

Hence, by (6.2), we have 
l open v adr Dato) Rx, SÑ 
ted Rjw de vw zę AREA Ga 
D"R : 

: „U 6.15 
(R, w (1 — R,w) ( ) 


wherein ' denotes derivation with respect to s. 
R,= X,— z, (s) 
We determine sy from the equation 
RR=0 KR >0 (6.16) 
With the notation 


Rw =n (6.17) 


we can rewrite (6.15) in simpler form: 


Ar DERNE 
4" us E | : ) | a (6.18) 


Eq. (6.18) was first derived by Mathisson (Mathisson 1942), by expanding the solutions of 
Eq. (6.1) in a series in powers of l/n. 


VII. Proof of the equivalence of the two methods 


We now proceed to prove that the formulas for the electromagnetie potentials of a 
dipole derived in Sections 4 and 5 are equivalent tot those derived in Section 6. We shall 
cousider the cases of the electric and magnetic dipoles separately. 
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a) Electric dipole. 


Let us investigate the form of the tensor D"' in the system wherein the dipole moves 
with velocity w. In the system at rest its form is that of Eq. (6.5). Introducing the notation 


Diss oss. pie. (7.1) 
we have, from the transformation properties (or from Eq. (6.6)): 
(D pep p (23) 
With the foregoing expressions, we can rewrite Eq. (6.14) as follows: 
|.(Dxvyx R- RD (Dxv)x yv - yD 
yin ym : 


A 


(Dx v) x R+ RD 
yos 


(I+R-ay — ly*va) 


D -R D EA D-R 
e gra dm a a m ss | a2 R q — Iys. 
p 272 WIĘ | 373 (14 y R:a ly (w a)) (1.3) 


In order .o compare (7.3) and (4.7), (4.9), we have to find the relationship between the 
vectors d and D. On writing the two formulas for the rest system, we see, on inspection, 
that 

d, — D, (7.4), 


It is readily proved, by considering the transformation properties of either quantity, that 
in an arbitrary system, 


D — yd (7.5) 
Taking into account the foregoing result, we see that (7.3) is identical with (4.7) and (4.9). 
b) Magnetic dipole. 
We introduce the notation 
M = (D3?, D13, D21) (7.6) 
Clearly, now 
(DH DŹ DR) = —Mxv (7.7) 
Eq. (6.14) now assumes the form 


MxR — R(Mx5) , Mx R— R(Mx v) 
~— s y3l3 


A= (14-5? R : a— bo - a) 


(MxvyR (Mxv)R 
A ES y? yam 


(l y? R a — Ly*v - G) (7.8) 


As above, in order to compare (5.19) and (5.21) with (7.8), it is necessary to find the rela- 
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tionship between the vectors M and m. Considerations identical with those sub (a) yield 


the. relationship 


M — ym (7.9) 


Substitution of (7.9) in (7.8) and simple transformations yield immediately (5.19) and (5.20). 


VIII. Generalization of the formulas derived 


It will be remembered that, in the case of a unipole particle i.e. one whose electric and 
magnetic moments jointly vanish, the eurrent density can be written as follows 


pa ji u" (s) d(x“ — z"(s) ds (8.1) 


Substitution of (8.1) in Eq. (6.3) yields Eq. (2.15) (see, e.g. Thirring 1958). We once more 
invoke the expression of the current density for a dipole: 


j= Ji D*'(s)0, ó(a”- : z" (s) ds (8.2) 


Eqs. (8.1) and (8.2) present the possibility of generalization. 
We introduce the notion of 2"-pole by prescribing its current density 


1 
Rx E. wo... 0,1 9 zy <Ą AL si 
juo pacc [EDS si Do, +++ Ba,9„0(x*— z"(s))ds (8.3) 


Quite generally, the expression for the current density of the singularities under considera- 


tion will have the form of the following series: 


= 2m (8.4) 


m/o 
Let us now revert to (8.3) in order to investigate the form of the tensor Di ?771 in more 
detail. First of all, we have to postulate that the current (8.3) shall fulfill the continuity 
equation. This leads to the condition 


HvGy... C vHO,.. € 


D m 1 D I m1 (8.5) 


Clearly, the index v is but apparently privileged, as any of the remaining indices could take 


ils place. We concentrate on the index » in order to fix ideas. Now, from the form of the 


subintegral expression, we see that the tensor D 7 ^7-1 can be assumed to be symmetrical 
in the indices o... ... gt i 
p n: 7m—1) 3 proc Im-ı (8.6) 
Moreover, let us note that the additional condition 
DW Ug = D e Us =... = 0 (8.7) 


can be prescribed. If the latter were not fulfilled, then, by applying the equality 


gr d 
u ar | (8.8) 
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and integration per partes, all terms not fulfilling it could be reduced to terms of 
lower m. 

We now proceed to determine formulas for the electromagnetic potentials of the 2” 
jeld. By substituting (8.3) in Eq. (6.3) and integrating m times per partes, we obtain in 
nalogy to Section 6 


— ])m wo. 6 
Ar = N" | pracom, ...0,0,Dretat— ze(ods 89 


+ 


By placing the differentiation operator before the integral and utilizing the function Dre 
in the form (6.2), we can obtain Eqs. (8.9) in the form 


rer ZRT: Be (8.10) 


m! n 


6 


In order to obtain insight into the physical meaning of the tensor D^ "m-1_ let us write 


Eqs. (8.10) for the multipole when moving and at rest. We obtain 


AF = 1a ES 66 Ok, OR 3 Dk... (8.11) 


whence we see immediately that 


ukky...k ME 


D É (8.12) 
is the electric moment of the 2" — pole, and 
pma (8.13) 


the 2" — pole magnetic moment (see, e.g. Stratton (1941)). 


IX. Relation to multipole radiation 


In order to exemplify the use to be made of the formulas derived here, we shall compute 
the potentials of a 2"-pole at rest, within the wave zone, i.e. taking into account only terms 
vanishing at infinity as //R. From the formulas 


Rae u (9.1) 


n=-ItRw=-1 (w=0) (9.2) 
we obtain, within the wave zone, 


(= Lye Ro, u Ro, Ry d" VO.. C. 
Am) = — Eur gp Dna (9.3) 


4 The Latin indices assume the values 1, 2, 3. 
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Taking into account (8.7), we have 


fi dg (— AM Ran, > Fa, Ry d" Dia Am + 
(m) m! Rma dt" 
(— 1)" Ra, Ru dm DM (9.4) 
m! Rm dt" 
(— 1)" Renat Ra Rk dn "o IR 9.5 
aT TEER d dt™ = i 0) 


Let us now separate these expressions into parts containing only the electrie or magnetic 


moment. respectively: 


(—1)" Rz SL dm 


z Ala) = a Im Pen; BPim) = Dem (9.6) 
° = R, c Rp, R; dm : 
MA = ( p = eH dm Dh m- M MO) —0 (9.7) 


These formulas can be considered to be generalizations of the well-known Rubinowicz- 
Blaton formulas (Rubinowicz and Blaton 1932). They yield the radiation of a time-variable 
electromagnetic multipole. Similar formulas have already been derived by various authors?; 
however, in each case a special from of the time-dependence of the eniire process was assum- 
ed, i.e. proportionality to e'”. The expressions of Eqs. (9.6) and (9.7) are free of this rest- 
riction. 
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DIELECTRIC SATURATION THEORY FOR MOLECULES WITH TWO 
AXES OF INTERNAL ROTATION 


By E. KLUK 
Institute of Experimental Physics, A. Mickiewicz University at Poznań. 
(Received April 27, 1961). 


The present communication is intended as the continuation of results obtained by A. Pie- 
kara, S. Kielich and A. Chełkowski in the field of dielectric saturation in rotational isomers and 
published in “Archives des sciences" (1959), and consists in a generalization of their concepts to 
cover rotational isomers presenting three carbon atoms in a normal carbon chain involving the 
existence of two axes of rotation. The essential part of the computations is carried out for the 
case of the CH, groups at either extremity of the chain containing a Cl, Br, or similar substituent 
to replace one of the hydrogen atoms, giving rise to a rotating group with a dipole moment. The 
present situation is more highly complicated than in dealing with a 1, 2-dihaloidoethane, inas- 
much as the number of rotational isomers is much greater and their energy states considerably 
closer to one another. Thus it seems reasonable to proceed by constructing a continuous potential 


barrier illustrating the interaction between the groups during rotation. 


In computing the dielectric saturation in compounds of the type of 1,3-dihaloidoprop- , 
ane, the formulas derived by A. D. Buckingham (1955) and A. Piekara and S. Kielich 
(1958) for gases are useful: 

Ku Kus) | pa 
= ao d BEES G l 
Ae ian | 30 43 T3 Eş d(e) (1) 
wherein (A. Piekara and 5. Kielich 1958 


IE, of F 
Ag Ë es ch 


and z 
| 9E, , 
e— ada | e DU A (2) 
with E  — denoting the external electric field strength, 
E, — the internal electric field strength, 
& — the strong field dielectric permittivity, 


(845) 
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u, -~ the momentary dipole moment of the molecule, 
$ 

a -its polarizability, 

N - the number of molecules per unit volume, 


V(t) — the intramolecular interaction potential, and 
dr —a volume element in configurational space. 


Thus. no specific intermolecular coupling is introduced, the only one to be retained 
being that occuring through the agency of the local field. This restriction seems justified, 
as the instability in the structure of the molecules under consideration is such as to affect 
intermolecular coupling adversely, so that the various kinds of intramolecular coupling’ play 
the decisive role. The entire problem thus reduces to one of finding the intramolecular 
potential and of computing the average dipole moment of degree 2n for the single molecule. 

The choice of the potential barrier adequately approximating the real interactions is 
an involved matter, as the otherwise universally applied method consisting in expanding 
e.g. cosines and sines of the angles determining the positions of the dipoles implies a too 
great number of constants, even if none but the initial terms are taken into account. Hence, 
it is necessary to assume that the interaction between the rotating groups is of dipole char- 
acter: 


V = kTY (cos ©, y — 3 cos ©, cos 6%) (3) 


wherein Y describes the energy of this interaction in kT units, 
Ois the angle between the dipoles, and 
5 


Oi, ©, -— those between either and the straight line connecting them. 
With these simplifications, the computation of the mean value < uż”) becomes somewhat 


less involved. On introducing the coordinate systems (a, y, z), (a, y" z), (a. y", 29 


in such a manner that the z-axis shall halve the angle 9 subtended by the two axes of rota- 


n A2) 


Fig. | 


tion, that z” and z” shall coincide with the first and second of these axes respectively, and 
that the planes (x, z), (x, z^), (a, z”) shall coincide with the plane of the carbon chain 


(Fig. 1), and using matrix notation, the following expressions are obtained for the dipole 
moments of either group: 


sin Ó cos p sin Ó cos y 


M, = m, | sin ô sin p]; Mę=m,|]sinó sin y 
cos Ó cos Ó 
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Herein, 44, and mp, are the absolute values of the dipole moments of either group, ó — the 
angle between axis of rotation and dipole moment, and p, V — the azimuth of each dipole 
in its respective coordinate system. 

The transformation from the system (w, y’, z^) and (a, y”, z”) to (x. y, z) are of the form 


D O (2, Ü 

| cos 5 W = s > COS > O sin > 
eo 1 0 

pz „GK 9 Ü 

| sin 5 U eos 5 | - sin 5 0 cosg 


By simple computations, the potential is readily expressed by means of the variable angles 


9 and y: 
> : mua Site. M 
= pi |= | 1 +00s? 5 sin? Ó cos p cos y + sin? Ó sin g sin y 
- 1 sin 9 sin 26 (cos p — cos | (4) 


Similarly, the square of the momentary resultant dipole moment is obtained: 


REM, [1+ 8? 28 cos 9 cos? ó + 2B sin? ô (cos 0 cos @ cos y+ sin @ sin y) 4 


+ B sin 9 sin 20 (cos p — cos y)]. (5) 
with 8 = ron averageing uż” over the angles p and y, we have 
Mis 
2n 27% £ V( y) 
f f o. y)e RT dg dy 
2n 0 0 
s s = Zn Za VU Fri š (6) 


ibis eL dą dy 
0 0 


On substituting herein the explicite expression for ję (p, y) and V (p, y), it will be seen 
that the computation of the foregoing integrals present considerable difficulties. These are 


overcome by using the transformation 

WERE Mew (7) 
The new variables are contained within the intervals 0 < £ < 4c and — 2z < £ < 2a, or, 
with respect to periodicity of the trigonometric functions, 0 < Ë < 4a. Since this transform- 


ation fails to yield directly the sines or cosines of the entire angles £ and £ in place of the 


expression cos p —cos y, it is necessary to use the following approximation: 


cos p — cos y = + (1 — cos £) (1 — cos I)? = 


il 
=. | = Zoos 5 =-= 2 cos à — TĘ cos? 6 + E Gos S cos ¿ — 16 cos? Ć + = (8) 
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to be taken positive for 0 < $ < Qn and 0 < € < 2z or 2z < £ < 4a and Za < Ẹ < 4n, and 
otherwise negative. This leads to iene analytie representations of the subintegral 


function in the various sectors of the interval of integration, and 


4141 an an 

"n g(cos £. cos 6) déd = 8 N gı(cos é, cos Ë) dë dË + Il ga (cos &, cos pasaż | 

ò 0 0 o 0 o 
Here, obviously, g,(cosé, cos Z) is the analytic form £ the function g (cos £, cosć) in the 
interval 0 < £ < 2m, 0 < ¿£ < 2z or 2z < £ < Az, 27 < C£ < 4r, whereas gą(cos &, cos C) is 
the one in O<&<2n,2n < ç <4a or Za LE< f 0 < £ < 2a. On substituting (7) 
and (8) in eqs. (4) and (5) and Me only the terms linear in the cosines of the angles 


£ and € we have 


Vio 
DAT A: 2 | 71,2 ros & + MA once 2 = a "OS ⁄ 08 C 
rA +Y,+!} P? cos $ ) 12 cos Ć (Wars [1 + b,,ə cos Š + c4, cos 6] 
wherein 
LES M 
Y, = re Y sin 26 sin 2 
ë J 1 
yi2 = } 1+ 1 cos? —| sin? ó F— sin 26 sin 0 
z > 8 
1 Ü l 1 
) 1,2 = Yi — C TÉ — Si 2 + — s = 
: 9 e 9 sin? ó 4 8 in 20 sin ? 
aig = m? [1 + p? + 28 cos 0 cos2ó + B sin 9 sin 20] (9) 
N 2 R dos ums Po oW E. 
^ia" 4,3 = Mis 2p sin“ ó sin*— $ p sin 9 sin 20 
3 ES T. 
C12 ` U2 = Mis | 2P sin? 6 cos? — F 4 B sin 9 sin 26 


With the notation already introduced, expression (6) is obtained in a form accessible to 


direct evaluation: 


; E HN ) tet te^ B, (YDB (YD + alby cze B, (Y3)B(Y3)] 
< 1o Ts s=0 
> ew BY DB, (YD + e B(Y3) B, (Y) 


(10) 


with 
Bim ji cos” ne” cos "dy 


0 
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In order to facilitate the discussion thereof, it is useful to simplify further the notation of 
eqs. (1) and (2) by introducing the correlation factors R, and R, as follows: 


2 2 4 
SED KW) — BC m$) 
Ręce e Rex TY (11) 


2 
HsM 


gut : : : . 
cis" configuration, being the maximum dipole moment of 


Herein, m= |t Heo] in the 
the molecule. Thus, eqs. (1) and (2) assume the form 


: DE, 
as l= daN E * E | FAL" (12) 
Now, on chosing one of the several known models of the local field, the hitherto undeter- 
mined parameter Y can be fixed form the polarizability and dipole moments of the isolated 
molecule. On the other hand, the first of eqs. (12) yields the dielectric saturation directly 
for various special cases. If the formulas (10) and (11) derived here are simplified so as to 
make the two axes of rotation coincide, forming one axis in common for both groups. 
Mizushima's factor R, for 1,2-dihaloidoethane 


TY) 


R,-1- AO (11a) 
and the factor R, derived by Piekara, Kielich and Chelkowski (1959) 
1 a) 5 [ B La) 
alle. on 


result immediately. Herein, /,(Y), L(Y) and 1,(Y) are modified Bessel functions. The 
relationship expressed by eqs. (lla) and (11b) for the case of 1,2-dihaloidoethane can be 


rendered still more exact if the potential barrier 
V — kT[Y' cos  4- Y" cos 39], 


better adapted to reality is introduced. The quantities of eq. (6) can now be evaluated if 


Da (= A Gs (Y) LY") 


(cos np > = = (13) 
> (1 Goro I; (Y) LY") 
pst 
has been computed, with 
- 1 1 
n ) | = IC TI Oi Neve hee FRE ss 
Gorn x ölp+ ( I) 3r+ ( 1) n] Tre Oop P: du 


Qs B=0 
As a matter of fact, the series of eq. (13) are sufficiently convergent for the quantities Y 


and Y" as occurring in practice. 
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Although the entire discussion hitherto was for a gas, it may be reasonably expected, 
when applied to liquids, to yield a dielectric saturation compatible with experiment. Indeed, 
on assuming the Onsager model in eqs. (12), the figures assembled in Table I are obtained. 
The measurements and computations are for room temperature. The Lorenz model, on the 
other hand, leads to entirely unrealistic results, as the value for KR, is then negative. It will 
be noted that there is agreement in sign between the experimental and theoretical data 


for dielectric saturation. 


TABLE I 
E = s Saar | os ra | 4 
| | | =%4l8|8 1&8 
| | | m= Me kTY-104|% a | w 
Gompound | N. 10-22 7 I3 = Hs E | © ea 8 e 
in D erg Rz „I B I 
| | KIER 
la. IN = = 
EZ — u u ps > | | — 
1,2-dichloroethane | 075 | 1.44 | 18 | 106 3.6 112,0 --17* 
1,2-dibromoethane | 0.69 1.54 | 1.80 4.9 9.7 +78 | --14 
| = - 
1,3-dichloropropane 0.62 | 1.44 | 2.04 | 9.8 5.6 LU —20 
1,3-dibromopropane | 0.59 1.52 | 2.04 9.2 10.1 EI IO 
1.4-dibromobutane | 0.50 1.50 | 1.80 7.8 | 3.6 +18.0 +2 


Maybe the too low theoretical values (except for 1,4-dibromobutane) could be brought into 
better agreement with experiment by altering the potentials assumed. Besides, the following 


ratios 

| g ÜHEOR | ¢THEOR 

A ECH;,CI, = A EC,H,Br. 3 
- = 0.71 — 0.56 
| EXP | „EXP 

A ECH.CI A SC,H.Br, 

| e THEOR THEOR 

C ECHC _ gą] CfGHB, o 41 

ec nci, A EGH.Br. 


go a long way to suggest that the discrepancies are due i. a. to the model of the local field 
being vastly inexact. It should be stressed that the latter plays a first rate part in the evalu- 
ation of the dielectric saturation. 

In carrying out the computations for 1,4-dibromobutane, the carbon chain was dealt 
with as being permanently in the “trans” configuration, whence the saliently high value 
of the saturation. As a matter of fact, however, liquid 1,4-dibromobutane contains a consid- 
erable percentage of “gauche” type configurations resembling that of 1,3-dihaloidopropane 
and providing a negative contribution to the resultant saturation. It is precisely this phe- 
nomenon that accounts adequately for the experimental decrease in saturation as the carbon 
chain grows in length, starting. from 1,3-dihaloidopropane. 

The author wishes to thank Professor A. Piekara, Director of the Institute of Exper- 
imental Physics of the A. Miekiewicz University at Poznan, for directing the present inves- 
ligation and for his valuable help. 
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APPENDIX 
Evaluation of the integral B. (Y). 
B,(Y) =| cos” ger ces 9 do = c — = Ip) 
; ayn 


The following relationship holds for the modified Bessel functions: 
: Pees, 
T, GY) =Ih-1(Y) = y ho) 


E (Y) — dn i (Y) = EY) 


5) 


whence 
0 


zy (7) = er, 1 (O| 


With the latter relationship, it can be proved that 


2m = 
= 7 
Jym ym In( Im > NIE n—m+2s ( Y ) 


whence 


Ba (y =s =. 2m =>, (”) lese) | 


s=0 


At higher values of the index m, the recurrential formula 


m! 


Bm+a(Y) = B,(Y)-4- (-Y)- 9*9 S'( Pr YP Bpsa (¥) 
p=0 


should be applied. 
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ON VARIATIONS IN THE DOMAIN STRUCTURE AND DIELECTRIC 
PERMITTIVITY OF BaTiO, SINGLE CRYSTALS DUE TO 
ONE-DIMENSIONAL PRESSURE 


By J. PIETRZAK 


Institute of Experimental Physics, A. Mickiewicz University, Poznań. 
(Received April 27, 1961) 


Polarized single BaTiO}. crystals were repeatedly subjected to one-dimensional pressure of 
200 kG/cm? applied along the [001] direction at room temperature. Irreversible changes of the 
domain structure and dielectric permittivity were observed. The single crystals ultimately pre- 
sented laminar domain structure of high stability and a permittivity no longer depending on the 


pressure. 


Introduction 


Inhomogeneous pressure has long been known to give rise to changes in the domain ' 
structure of barium titanate (Blattner, Känzing and Merz 1949, Forsbergh 1949); however, 
the trend of these changes in the presence of considerable, homogeneous mechanical stress 
has hitherto failed to be clarified. 

Ksendzov and Rotenberg (1959) found that the variations in dielectric permittivity of 
BaTiO, single crystal due to one-dimensional pressure depend on the direction along which 
it is applied. Stankowski (1960), working in this laboratory, found the changes in domain 
structure of BaTiO, single crystals produced by one-dimensional pressures of up to 10 at 
to be reversible. On the other hand, the present author found that higher pressure gives rise 
to irreversible changes in the domain structure and dielectric permittivity. Microscopic 
observation of the domains was carried out on removing the single crystal from the press. 

It was the aim of the present investigation to clarify the effect of one-dimensional 
pressure on the domain structure and dielectric permittivity of barium titanate c-domain 


single crystals. 


Experimental 


BaTiO, single crystals grown by Remeika'a method, in the shape of thin plane-parallel 
plates 0.1—0.3 mm thick, were used. The crystals were polarized in an electric field of 
about 5 kV/cm applied in the direction [001]. Monocrystals for optical investigation were 
polarized with liquid electrodes, whereas those intended for measurements of the permit- 
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tivity were provided with silver electrodes by evaporation in vacuo. One-dimensional pressure 
swas applied in a press. Observation of the domain structure was carried out with a polari- 
zation microscope. The capacitance of the single crystals was measured with a Schering 
bridge at a frequency of 2 Ke/sec and a measuring field of about 10 V/em. Observation of 


domain structure and dielectric measurements were effected at room temperature. 


Results 


Fig. | shows the variations in the domain structure of a c-lomain single crystal arising 
from one-dimensional pressure applied in the direction [001]. The position of the crystal 
on the microscope support was always such that one of the cubic axes subtended an angle 
of 45° with those of the crossed polarizers. Observation of the domain structure was carried 
out during the 15 minute intervals between applying one-dimensional 200 kG/em? pressure 
for periods of 15 minutes. In c-domain single erystals, one-dimensional pressure gives to 


domains whose direction of polarisation is perpendicular to that of the pressure applied. 


Fig. 1. Domain str i i j 
ea] e b uc e x rle ervst; cte ^ inr S 1 Š 
sę wa structure in a single crystal aJ not subjected to pressure, b] subsequent to applying pressure of 
2 g/cm? for 15 mi cw E E x iis ‘ ea : 
= £ i ^» min., once, c) subsequent to applying pressure of 200 kg/cm? for 15 min., at intervals of 
o min., four times, d] subseque í ing 200 kg/cm? for 15 mi i "15 mi [ i 
d] subsequent to applying 200 kg/em? for 15 min., at intervals of 15 min., four times; ob- 


servation in monochromatic light. 
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After each period of pressure, the width of the domains is seen to diminish, whereas their 
number augments. A domain structure arises, presenting the form of narrow stripes differ- 
ing in width and running parallel to one of the cubic edges of the monoerystal. The polar 
axes of the domains (which will be termed laminar, to denote their structure) are directed 
with respect to one another as shown schematically in Fig. 2. In observing the monocrystals 
when one-dimensional pressure varied, it was found that propagation of thin, wedgelike 
laminar domains throughout the bulk of the crystal constitutes the only process by which 


ALAL a 


Fig. 2. Laminar domains and directions of polar axes (—]) in single crystal on applying one-dimensional pressure 
several times. 


new domains arise. 


The system of such domains exhibits a high degree of stability. No noticeable changes 
were observed to occur in the configuration of laminar domains during several weeks. Neither 
were any permanent changes found to result after an electric field had been applied, nor from 
heating the monocrystal several degrees above the Curie point and then cooling it to room 
temperature. However, if the crystal is heated at 700^C and cooled to room temperature, 
the configuration of the laminar domains changes, the crystal remaining multi-domain. 

~ Fig. 3 shows the relative variations of the dielectric permittivity in a Ba TiO, single 
crystal versus one-dimensional pressure applied parallel to the [001] direction and to that 
of the measuring field. Graph a represents = versus the pressure for a “virgin” single 
crystal. The result is in agreement with measurements by Ksendzov and Rotenberg (1959). As 
the crystal is subjected to the one-dimensional pressure for the second, third, etc. time 
(Fig. 3, graphs b, c, d, e), the consecutive relative variations of the permittivity are found 
to become progressively lesser, and finally (Fig. 3, graphs d, e) altogether ceäse to occur. 
The foregoing variations are positive, corresponding to an increase in the permittivity. 


DE: coc reru 
The values of — in Fig. 3 are referred to the permittivity of a “virgin” single crystal. 
8 


Discussion 


The very considerable stability of the laminar domain configuration arising in a single 
crystal subjected to one-dimensional pressure is a clear indication that, m addition to the 
laminar domains, local deformations of the crystal are given rise to. This is corroborated by 
the fact that, on having been subjected several times to one-dimensional pressure, the 
monocrystals exhibit optical anisotropy due to internal stress in polarized monochromatic 


856 


100 200 


ees 


€ | n 3 
Fig. 3. — versus one-dimensional pressure (graphs a, D, c, d, e resulted from consecutive application of pressure 
E 


a number of times). 


light (Fig. 1 d). Therefore, due to pressure, the single crystal undergoes plastic deformation 
as a result permanent laminar domains structure (twinning) and dislocations with the slip 
planes lying along the 90°-walls. 

The presence of unusually stable laminar domain configuration is corroborated by 
measurements of the dielectric permittivity versus the presure (Fig. 3, graphs d, e), as, 
in the crystal deformed, the permittivity no longer depends on one-dimensional pressure. 
This proves that the latter no longer raises the number of domains. Thus, the monocrystals 
can be said to be tempered as a result of plastic deformation. 

The author wishes to thank Professor A. Piekara for his direction throughout the 
present investigation, and for his valuable hints and advice. 
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GROUP TRANSITIONS IN CRYSTALLINE PHOSPHORS 


By T. I. KAKUSHADZE 
A. 5. Pushkin Pedagogical Institute, Tbilisi, USSR 
(Received December 12, 1960; revised manuscript received May 6, 1961) 


The present paper attempts to apply the group transition theory to the problems of the 
luminescence of crystalline phosphors!. 


l. Group Transitions? 


In the Peierls conductivity theory, the wave vector of a single electron changes by 
a whole vector of the reciprocal lattice. Starting from the work of Bloch [7, 8], Houston 
showed that the Hamiltonian of the system electron-lattice remains invariant if the wave 
vector, not of a single electron as assumed by Peierls [9], but of a large group of electrons ` 
changes by a whole vector of the reciprocal lattice. For a particular case of pair transitions 
Houston proved that the matrix elements containing the products of the interaction of two 
electrons with the lattice vanish unless the following condition is satisfied 


(RR eK, Ke) een. (1) 


Here n is an integer, E K, and Ri ie, are the wave vectors of the electrons before and 
after the transition, and G is a whole vector of the reciprocal lattice. 

For the particular case when n=0, Eq. (1) is the usual pair transition (without Peierls’ 
umklapprozesse) 


KOR ZR eon (2) 


If eqs (1) or (2) are fulfilled, the matrix elements differ from zero and consequently 
the corresponding transition does occur. 

The energy band scheme of the semi-conductor makes it possible to give a graphic 
interpretation of a group transition with a subsequent excitation of an electron of the va- 


1 The application of the group transition theory to conductors has yielded fairly satisfactory results [1-5]. 
2 Ip a somewhat briefer form, this work is now being published in the “Annalen der Physik” [Berlin]. 
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lence band into the conduction band (Fig. 1) (1 and 2 in Figs. 1-4 designate the valence 
band and the conduction band respectively). The arrow aa' corresponds to the energy Eaa’ 
of the transition of an electron from the valence band to the conduction band. The arrow 
bb’ corresponds to the change of the energy of the system due to all other electrons in- 
volved in the group transition (Kaa' = Ebb’). 


The total number p, of group (pair) transitions of conduction electrons into the valence 
band (Fig. la) per time unit is proportional to the product of the number N of conduction 
electrons by the square of the number of holes in the valence band m? and the number M 
of the electrons in the lower part of the valence band (N, being independent of tempera- 
ture), $.e., 


Py = AN,Nm? = G,Nm?, (3) 


where G, is a constant. 

At the absolute zero the thermal energy of the system of the electron is equal to zero 
(m=N=0, Eq. (3)). Therefore at O°K only such group transitions occur in which the 
electrons of the ground band simply exchange places, without changing the physical state 
of the system. 

At other temperatures Fermi holes appear in the valence band. Every group transition 
leads to one sufficiently deep hole b” (Fig. 1b) due to a large number of Fermi holes [3]. 
If there is a hole b” then an energy conserving (isoenergetic) transition of an electron from 
the valence band into the conduction band (Fig. 1b) may occur. It is essential that both 
the formation of a deep hole and the transition of a valence band electron into the conduc- 
tion band occurs in a single process of a group transition. A large amount of energy released 
by the group of electrons is absorbed by one electron of this group, ensuring thereby the 
energy conserving transition of this electron from the valence band into the conduction 
band. Due to energy conservation in group transitions the thermal energy of the electrons 
remains unchanged. 

A strong perturbing influence due to group transitions led Houston to the conclusion 
that the group transitions must have higher probabilities than single transitions. In partic- 


ular, pair transitions occur much more frequently than single transitions. According to 
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Houston [6] the contribution, to the perturbation, of the term proportional to the first 
power of the displacement of the ions is very small. Because of the symmetry of the inter- 
action between an ion and an electron, the average of the first derivative will vanish. If the 
second derivative is also taken into consideration, quadratic terms in the displacement 
must be included. 

As a result of a group transition an electron of the valence band may also pass to the 
conduction band when the width of the forbidden gap is so large that the thermic excitation 
of the electron to the conduction band is rather improbable according to Fermi-Dirac 
quantum statistics. 

As was shown by Houston [6] group transitions arise due to strong electron-lattice 
interaction?. 

In group transitions, the law of conservation of wave vectors in the non-oscillating 
lattice, without taking account of Peierls’ umklapprozesse, is of the form 


> 


"n ay n 
> K, = > i? 
i=] emi 


e.g. for a pair of electrons: K,--K, — K, -- K5, where K, and K, are the initial, and Ki and K; 
are the final wave vectors. Hence the energy of the pair after the transition will be 


om . 2m 2m . 2m 2m 2m 


2 2 
: LIE RE eee [e cosa 2KIKš cos z) 


where & and «' are the angles between the wave vectors before and after the transition re- 
spectively. The energy of the pair is liable to increase or decrease considerably with the 
angles œ and g’. A group transition can be interpreted as consisting of a large number of ' 
such transitions. The energy conservation of a group transition requires that a certain number 
of transitons should give a posilive expression in the brackets and a certain number a ne- 
gative expression. The positive expression in the brackets corresponds to the production 
of a deep hole in the space of wave vectors (Fig. 2a) and the negative one the disappearance 
of the hole (Fig. 2b). Strong electron-lattice interaction leads to a large number of simulta- 


Fig. 2 


neous processes of both types. Such group transitions with conservation of energy lead 
in final analysis to the transition of at least one electron to a high vacancy level. 


3 An obvious experimental proof that group transitions have high probabilities is the experimental exci- 
tation of luminescence by electrons of small kinetic energies or by hot ions (see below). 
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Apart from group energy conserving transitions, impurity atom donors may be the 
carriers of electrons to the conduction band. Similarly, the origin of part of the holes in 
the valence band is due to the presence of impurity atom acceptors in the semi-conductos 
[12]. At low temperatures the number of conduction electrons and that of holes in the va- 
lence band depend on impurity atoms, whereas at sufficiently high temperatures the change 
of the number of holes in the valence band and that of electrons in the conduciion band 
are entirely due to inter-band transitions. 

It has been established [12] that the temperature change of the number of holes in the 
valence band can be described by the formula 


Eeff 
mx Ae UT > (4) 
where E,g is the effective energy value, equal to a few hundredth of eV (the energy of the 
acceptor level being counted from the top of the valence band) at low temperatures, and 
to several eV (the energy width of the forbidden gap at sufficiently high temperatures) at 
room temperature and above. 
According to (3) and (4) the mean time 7, of the group transition of a conduction 
electron into the valence band is equal to 
1 Eeff x 
a= > = pe 2kT 5 
Un > ; (5) 
where g is a new constant and Gym? the probability of transition of a conduction electron 
to the valence band. 


2. Transition Times Tq, Tp, Tr 


The time 7, of the non-radiative transition of a conduction electron strongly depends 
on the band structure of the semi-conductor. Thus, in a semi-conductor (Fig. le) with | 
the width of the valence band less than that of the energy gap there are no transitions of 
low multiplicity (double, triple transitions, ete.) which usually occur much more often 
than those of a higher multiplicity. Therefore, the time 7, of the transition of an electron 
from the conduction band into the valence band should be much larger in semi-conductors 
of the latter type than the respective time Tọ in those with a broad valence band (Fig. la, b). 

It can be supposed that the probability of a group transition strongly depends on the 
magnitude of the momentum p of the particle hitting the semi-conductor. The photon- 
induced transition of an electron is essentially restricted by the selection rules [12, 10]. 
A practically no-momentum photon (p = hy/c) is unable to effect the transition of an electron 
from one state into another in the same band if the final state is occupied; owing to its low 
momentum, the photon cannot cause a group transition of electrons. 

An entirely different effect should be expected if a semi-conductor is bombarded by 
hot electrons and ions of a gas in a flame into which the sample under test is placed. Having 
momenta 10* to 10% times larger respectively than the photons with the same energies, 


D; mwc 2c mw? 

ae ch TR" at h = ET where p; and py are the momenta of an electron or an ion and 
Df w v 

a photon respectively. 
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hot electrons and ions loosen the valence band. The mechanism is represented schemati- 
cally in Fig. 3. 

Let the incident particle interact with a valence electron a to which the energy Eaa’ 
may be transferred. If the state a’ is filled, a large momentum of the incident particle causes 
perturbation, leading to an energy-conserving group transition, shown by the arrows bb’ 
and cc', with the subsequent release of the state at the level a/, A large momentum of the 
incident particle ensures thereby the transition of an electron a to the a’ state (with the 


————————— 


2 


25 


resulting deep hole a). The probability of such a transition differs from zero (see below). 
We suppose that a hot ion may cause the production of a deep hole, since we admit that 
a high momentum particle is capable of producing a group transition. Consequently, in 
a group transition the aa’ transition may occur for the excitation by a photon even if the . 
latter is forbidden by the Pauli principle (the state a’ being filled). 

The bombardment of a semi-conductor by ions or electrons loosens the valence band 
completely. A group transition in such a conductor essentially increases the probability of 
the umklap of an electron from the valence to the conduction band (Fig. 1b). 

The laws of conservation of the sum of wave vectors and electron energies are fulfilled 
in the above processes because of the group character of the transition. 

The following considerations are essential in comparing the results of the group-transi- 


tion theory with experimental facts. 
At room temperatures Ty for different semi-conductors assumes the values[12-13] 


T OS ORF sec, (6) 
The rate of change of Tę with temperature is, according to Eq. (5), given by the formula 


dt, Be 
dT ~~ Ep 0 


The energy band sheme of an ideal crystal phosphor is known to contain (apart from 
the valence band and the conduction band) F and D localized levels situated in the energy 
gap (Figs. 4 and 6). The F-levels are known as the trapping centres of a conduction electron 
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and the D-levels as the luminescence centres [12, 14]. It is also known that a conduction 
electron, before passing to the luminescence level with emission of a luminescent quantum, 
can be captured repeatedly by the trapping centres with a resulting prolonged afterglow. 

The times Tp and Ty of the transition of a conduction electron to one of the deep localized 
Zzh$(R 
e? E°a? 
that the energies of the transition of a conduction electron to the D and F localized levels 
are equal to several eV and several hundredths of eV respectively, a~ 10-3, we find 


levels F or D are evaluated by the formula v — (dipole transition [10]). Assuming 


To~ 10-* sec,  Tg--10-? sec. (8) 


It should be emphasized that a group transition always occurs with energy conserved and 
is, consequently, radiationless, whereas the transitions of a conduction electron to the F 


Fig. 4 


and D localized levels are accompanied by emission of a photon. The transition of an electron 
will either occur with the radiation of a luminescence quantum or will be radiationless de- 
pending on whether an individual electron passes to the D local level (Fig. 3ab) or the tran- 
sition involves a large number of electrons (Fig. 3c, the transitions aa’, bb’ or cc’). 

Eqs. (6) and (8) lead to the condition 


T> TR Tp (9) 
or 

TpŻ ToS Tp (10) 
or 

TEŻ Tp> Tg. (11) 


The present paper is confined to the investigations of conditions (9), (11) involving 
group transitions of electrons. 

The following experimental fact directly illustrates the origin of group transitions as 
a result of the bombardment of a phosphor by microparticles with large momenta. Shrader 
and Kaisel [11] bombarded the phosphor ZnO [Zn] with a beam of electrons with kinetic 
energies not exceeding 0.1 eV, and the luminescence radiated contained photons with 
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energies from 1.8 to 3.1 eV. Such a production in a semi-conductor of a large quantum out 
of smaller ones is a graphic evidence in favour of the existence of group transitions. The 
incident electron interacting with the valence electron a (Fig. 3) must indeed excite it into 
‘the a’ state. Yet the release of the a’ state is necessary for this transition to take place, and 
this occurs as a result of a group transition due to a large momentum of the incident elec- 
tron. A group transition of this type is represented on Fig. 3 by the arrows bb’ and cc’. 

The production of deep holes similar to a may result in the group transition of a valence 
electron to the conduction band, as shown in Fig. 1b. The reverse transition of the electron 
to the D localized level yields the corresponding luminescence shown in Figs. 4a, b. The 
luminescence yield (ratio of the intensity of luminescence to that of energy absorption) 
in ref. [11] suggests that the probability of a group transition in the bombardment of a phos- 
phor by fast electrons is quite considerable. 

Another convincing proof that the bombardment of a phosphor by microparticles with 
large momenta is responsible for group transitions is the well-known experimental pheno- 
menon of supertemperature luminescence. The experiment reduces to the following. Certain 
semi-conductors (NbO, ZnO, ZnS, TiO, and others) are capable, when incandescent, of 
emanating supertemperature luminescence [15, 16]. Here the transition of an a electron 
to the a’ state with a simultaneous energy-conserving group transition (Fig. 3, transitions 
bb’ and cc’) is effected by hot electrons or incandescent ions with large momenta. If there 
are deep holes which have arisen in the above manner, group transitions with the ejection 
of a valence electron to the conduction band occur as shown in Fig. 1b. The transition of 
the latter electron to the D localized level is accompanied by due luminescence. It is the su- 
perimposition of the latter on the thermal radiation of a semi-conductor that leads to su- 
pertemperature luminescence, or candoluminescence as it is called, the intesity of supertem- 
perature radiation reaching large values in some case. The luminosity of niobium oxide 
in the blue region of the spectrum at 560°C is 85,000 times as strong as the respective in- 
tensity of black radiation [24]. The fact that the temperature of the sample is not very high 
shows that the kinetic energy of an overwhelming majority of hot ions is considerably less 
than a quantum of luminescence. This selective character of irradiation of certain phosphors 
warrants the conclusion that the probability of a group transition of electrons in the bom- 
bardment of the latter by hot electrons and ions is not small. These examples indicate 

that group energy-conserving transitions in solids may be caused by high-momentum 
electrons. The interactions of electrons with crystal-lattice ions may lead to the appearance 
of high-momentum electrons in the valence band. The origin of such transitions can therefore 
be expected as occurring thanks to internal electrons, without external bombardment [6]. 


Comparison of Theoretical Results with Experimental Facts 


1. “Anomalous” Contraction of % 


An essential conclusion of the group transition theory can be presented as follows. 
It was pointed out above that at room temperatures the time 7, of the group transition of 
a conduction electron into the valence band for different semiconductors assumes values 
from 10-1 to 10-12 sec (Eq. 6). According to (7) Tọ changes especially rapidly in the tempera- 
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ture range of 300°-400°K, when the donors and acceptors are practically empty and Ej 

becomes —2eV. An estimate made according to eq. (5) shows that the change in the BĘ 
perature of a sample in the above range corresponds to the decrease of the time 7, by a factor 
of 107. According to (6) and (7), therefore, zç for many semi-conductors must, at not too 
high temperatures, become less than the time necessary for an electron to skip one lattice 
constant. Contradicting the modern concept of the mechanism of the electric conductivity 
of conductors and semi-conductors, this fact has nevertheless been experimentally confirmed 
in full accordance with the group transition theory. In a number of semi-conductors, 
especially in ionic crystals, the electron free path is much less than the lattice parameter [12] 
Such a considerable decrease of the electron free path may be attributed to the participation 


of a large numbér of electrons in each group transition. 


2. Radiationless Transition 


'The group transition theory leads to another, still more important conclusion. Group 
transitions, as was pointed above, occur with the conservation of energy and without radia- 
tion. A large amount of energy released from the transition of a conduction electron to the 
valence band (Fig. la, transition cc’) is absorbed by another electron (or a group of elec- 
trons) of the valence band (Fig. la, transition dd"), which results in a deep hole d. As a result 
of the subsequent interaction of electrons with the oscillations of the lattice, the hole d 
reaches the top of the valence band, exchanging the large energy of the transition for 
a large number of phonons. That is the mechanism of the radiationless electron transition 
in the framework of the group transition theory. 


3. Luminescence of Crystalline Phosphors 


Inequality (9) is the condition that in the preliminary excited phosphor the transitions 
of conduction electrons to the F and D localized levels occur with the emission of a photon. 
Indeed, according to (9), before a radiationless transition of a conduction electron to the 
valence band can take place, its transition to the F or D localized levels occurs with emission 
of a photon (Fig. 3a). 

For a certain number of phosphors inequality (10) must hold. This inequality is the 
condition that the preliminary excited phosphor must glow. Indeed, according to (10), 
before the group radiationless transition of a conduction electron takes place, there may 
occur its transition to the activator level with radiation luminescence. The capture of the 
electron by the F-level, on the other hand, must according to (10), be radiationless (Fig. 3b). 
This luminescence may be actually observed in certain phosphors (ZnO, ZnS, CaS etc.). 
In particular, the substances ZnS, ZnO, CaS etc. begin to glow after being irradiated by 
X-rays or cathode rays. 


4. Candoluminescence 


If condition (9) or (10) holds up to high temperatures (the forbidden gap being larger 
than the valence band) bright luminescence must be observed corresponding to the transi- 
tions of conduction electrons to the D localized levels of luminescence (Figs 4a and b). 
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The vigorous penetration of holes into the depth of the valence band (Fig. 2) (necessary for 
the transition of an electron of the valence band to the conduction band) is ensured by the 
bombardment of the phosphor by hot electrons and ions in the flame in which the incan- 


‘descent semi-conductor has been placed. 


he Cola lare" 


The preliminary excited phosphor for which the conditions (11) hold at room tempera- 
ture cannot glow. True enough, before the transition, with emission of a photon, of a con- 
duction electron to the F or D localized level takes place, there occurs its group radiation- 
less transition into the valence band or to localized levels (Fig. 4c). 

As the temperature decreases, inequality (11), according to (5), must go over into the 
condition 


TED T0) Tp (12) 


and consequently the preliminary excited phosphor must begin to glow. According to (14), 
there occurs indeed, a preliminary transition of a conduction electron to the activator level 
with radiation of a luminescence quantum (Fig. 4b, transition b). This phenomenon — the 
irradiation of a phosphor at low temperature while it remains dark at room temperature — 


has been actually observed [14] and is known as “cold flare”. 


Adirovitch [18] observed a cold fłare in preliminarily excited and then irradiated sam- 
ples. When the sample was cooled to 80 °C, it “flared”, whereas it remained dark at room 
temperature. After the flare ceased the temperature of the sample was lowered to room 
temperature and then the sample was again cooled to the temperature of liquid air. Adiro- 
vitch observed from 3 to 5 flares on end without additional excitation of the phosphor. 
According to Adirovitch, 3, 4 or 5 cold flares lead to the emptying of small localized levels. 

As a result of prolonged rest, the sample became again capable of a cold flare. For exam- 
ple, 2 or 3 cold flares were observed a day later. It can be assumed that during the rest period 
a number of electrons passes from deep levels to shallow adhesion levels. In other words, the 
“cold flare” can, with allowance for deep and shallow levels, be interpreted in the frame- 


work of the group transition theory. 


6. Quenching of Luminescence 


As the temperature rises, condition (10) goes over, according to (10), into the following 


one 


te tp ty. (13) 


Under condition (13) the transitions to the F and D levels are radiationless. Before the 
transition, with the radiation of luminescence, of a conduction electron to the F and D 
localized levels can, indeed, take place, there will occur, according to (13), its radiationless 
transition to the valence band (Fig. 4c, transitions aa’, bb’ or cc’). This phenomenon — 
temperature quenching of preliminary excited phosphors — is actually observed, in full 
accordance with the group transition theory, and is known as the “quenching of lumines- 
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cence”. The total quenching of many phosphors occurs at a temperature ~400°K at which, 
according to (5) and (6), 79 must be much less than the time vp —10-* sec of the radiative 
transition of an electron (we assume E,g —2 eV). 

Ch. Peyrou [19] observed temperature quenching in samples of ZnS and ZnS.CdS 
(50% of ZnS and 50% of CdS) (Fig. 5a) and Garlick and Gibson [20] in samples of ZnS.Zn 
(Fig. 5b). In all these cases complete quenching of luminescence was observed at tempera- 


tures —400°K, in accordance with Eqs. (5) and (8). 


ZnS - Zn 
1 A*253.7 mu 
yt 2 A=313.2 mp 
1 3 A= 365.3 mu 
2 < A= 404.3mu 


93 133 173 203 253 293 333 373 413 TOK 100 200 300 400 500 T°K 


a b 
Fig. 5 


At a sufficiently low temperature condition (13), according to (5), goes over to the 
new condition 


TR 192» Tp. (14) 


Under condition (14) an electron in a preliminary excited phosphor passes to the D level 
with the radiation of a luminescence quantum if the thermal release of the corresponding 
localized levels is still possible at a given temperature. This phenomenon has likewise been 
observed in full accordance with the group-transition theory. 

It will be noted that inequality (10) is the condition of the “cold flare" of either usual 
luminescence or candoluminescence, depending on whether this condition is fulfilled at 


a low temperature (below room temperature), at room temperature, or at higher practically 
available temperatures. 
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7. Temperature Limits of Luminescence Bands of Multiactivator 
Phosphors 


A matter of considerable interest is the existence of the temperature limits of the lumi- 
nescence bands of multiactivator phosphors. Fig. 6 represents the band scheme of a two- 
activator phosphor. The recombination times of a conduction electron with the activator 
levels D, and D, will be designated by tp, and Tp, respectively. Then at temperatures for 
which the following inequality obviously holds 


T9 S Tp, ŚTp,: (15) 


in the luminescence spectrum of a preliminarily excited twoactivator phosphor, there will 
be no band corresponding either to the transition of an electron to the D, level or to its 
transition to the D, level. As the temperature decreases (at T'—7)) the following condition 


will, according to (5), hold 
To — Tp, ŚTp,: (16) 


Consequently, there will arise in the luminescence spectrum a band corresponding to the 
transition of a conduction electron to the activator D,-level. On the other hand, the band 
which may be attributed to the transition to the D, activator will be absent. As the tempera- 
ture falls still lower and the following condition becomes valid 


Tp, Tp, To (17) 


the second band emerges as well. The disappearance of either of these bands is connected 
with the energy of thermal release of an electron from the activator with subsequent tran- 
sition into the valence band. The first to disappear will be the band which takes more work 
of exit of an electron from the respective localized level. Both bands disappear at a suffi- 
ciently low temperature. 

Let in particular T, and T, be respectively the temperatures of the beginning of the 
thermal release of electrons from the D, and D, activators? and T, and T, be the tempera- 
tures at which 79 becomes of the order of vp, and Tp, respectively. Then the inequality 


Teste Tee Le (18) 


is the condition that in a preliminarily excited phosphor the first band will appear before 
the second, and both bands will exist in the temperature intervals T,—T, and T,— T; 
respectively. For different phosphors the quantities Ti, T5, T, and T} may exist in different 
ratios, which will affect the character of the origin of the temperature bands. In particular, if 


Tyś d e Ii (19) 


the bands will not be observed at all. 
This phenomenon — the origin. of luminescence bands with different temperature 


limits in the spectra of multiactivator phosphors — has been actually observed. 


5 T, and T, are determined by the work necessary for the release of an electron from the respective acti- 


vator and do not at all depend on the energy position of the local levels. 
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Different temperature behaviour in different phosphorescence bands in the same lumino- 
phoric compounds have attracted attention of many scientists. Thus, Schon [21], Kla- 
sens [22] and Adirovitch [23] observed the temperature behaviour of luminescence bands 
of the ZnS. Cu phosphor. At comparatively low temperatures they observed the band cor- 
responding to the transition of electrons to the D level of Zn, yet with rising temperature 
the Zn band disappeared and the Cu band grew in intensity. 

Thus the group-transition theory explains, at least qualitatively, the complicated pro- 
cesses occurring in crystalline phosphors (the excitation of luminescence by low energy 
electrons, anomalous contraction of a free path in certain semiconductors, radiationless 
transitions, candoluminescence, cold flare, quenching of luminescence, and the temperature 
limits of the luminescence bands of multiactivator phosphors). We believe that the further 
development of the group-transition theory will make it possible to study in qualitative 


terms other complex processes occurring in solids. 
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REVIEWS OF BOOKS 


ERWIN SCHRÓDINGER 


GEIST UND MATERIE 
Friedrich Vieweg und Sohn, Braunschweig 1956, 78 S. 


Dieses Buch enthält die Übersetzung der Vorträge, die Schródinger am Trinity College 

in Cambridge im Oktober 1956 gehalten hat, und die in der englischen Sprache unter dem 
Titel „Mind and Matter* publiziert wurden. In diesem Buche befasst sich Schródinger in 
der anschaulichen Form mit den Problemen des Verhältnisses von Materie und Bewusstsein 
und der lebenden und toten Substanz. 
Im einleitenden Kapitel wird das Problem der physikalischen Grundlage des Bewustseins 
| gestellt. Dann folgen die Betrachtungen über die allgemeinen Züge der Evolution des Le- 
bens und der Menschheit mit der Kritik der heutigen industriellen Zivilisation. Weiter wird 
eine der grössten Schwierigkeiten der heutigen Wissenschaft besprochen: das Verhältniss 
vom Objektiven und Subjektiven und das wissenschaftliche Weltbild, wo das „ich“ des For- 
schers keinen Platz hat. Dann bespricht Schródinger das Grenzgebiet zwischen physiologi- 
schen und psychologischen Erscheinungen; im Zusammenhang damit wird die Evolution 
des physikalischen und philosophischen Begriffes der Zeit vom Altertum bis zur Gegenwart 
geschildert. Im letzten Kapitel, welches den Titel „Das Geheimniss der Sinnesqualitäten“ 
trigt, kehrt Schródinger wieder zum Problem des Objektiven und Subjektiven zurück, 
indem er das Verhältniss der Messungsvorrichtung und des Experimentators darstellt und 
die Rolleeder Experimente in der Bildung der physikalischen Theorien schildert. 

Das Buch enthält viele interessanten Betrachtungen und Bemerkungen des grossen 
Physikers über die Probleme, welche eine grosse Rolle in der Bildung der Weltanschauung 
spielen und welche an der Grenze von Naturwissenschaft und Philosophie stehen. 
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